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ABSTRACT 

The type la supernova (SN la) rate, when compared to the cosmic star formation 
history (SFH), can be used to derive the delay-time distribution (DTD. the hypothet- 
ical SN la rate vs. time following a brief burst of star formation) of SNe la, which 
can distinguish among progenitor models. We present the results of a SN survey in 
the Subaru Deep Field (SDF). Over a period of three years, we have observed the 
SDF on four independent epochs with Suprime-Cam on the Subaru 8.2-m telescope, 
with two nights of exposure per epoch, in the R, i! , and £ bands. We have discovered 
150 SNe out to redshift z 2. Using 11 photometric bands from the observer-frame 
far-ultraviolet to thc^ near-infrared, we derive photometric redshifts for the SN host 
galaxies (for 24 we also have spectroscopic redshifts). This information is combined 
with the SN photometry to determine the type and redshift distribution of the SN 
sample. Our final sample includes 28 SNe la in the range 1.0 < z < 1.5 and 10 in the 
range 1.5 < 2: < 2.0. As our survey is largely insensitive to core-collapse SNe (CC SNe) 
at z > 1, most of the events found in this range are likely SNe la. Our SN la rate mea- 
surements are consistent with those derived from the Hubble Space Telescope (HST) 
GOODS sample, but the overall uncertainty of our 1.5 < 2: < 2.0 measurement is a 
factor of 2 smaller, of 35-50 per cent. Based on this sample, we find that the SN la 
rate evolution levels off at 1.0 < z < 2.0, but shows no sign of declining. Combining 
our SN la rate measurements and those from the literature, and comparing to a wide 
range of possible SFHs, the best-fitting DTD (with a reduced = 0-7) is a power 
law of the form ^{t) ex. , with index /3 = —1.1 ±0.1 (statistical) ±0.17 (systematic). 
This result is consistent with other recent DTD measurements at various redshifts and 
environments, and is in agreement with a generic prediction of the double-degenerate 
progenitor scenario for SNe la. Most single-degenerate models predict different DTDs. 
By combining the contribution from CC SNe, based on the wide range of SFHs, with 
that from SNe la, calculated with the best-fitting DTD, we predict that the mean 
present-day cosmic iron abrmdance is in the range Zpc = (0.09-0.37) Zpe,®- We fur- 
ther predict that the high- 2; SN searches now beginning with iySTwill discover 2-11 
SNe la at > 2. 

Key words: surveys - supernovae: general - cosmology: miscellaneous - cosmology: 
observations 
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1 INTRODUCTION 

Supernovae (SNe) play important roles in a variety of astro- 
physical settings, from galaxy evolution to the metal enrich- 
ment of the interstellar medium, as catalysts of star forma- 
tion, and as distance indicators. SNe are separated into two 
main physical classes: core-collapse SNe (CC SNe), which 
include all type II SNe (i.e., those objects which exhibit 
obvious H lines in their spectra) and Type Ib/c SNe (i.e., 
spectra lacking H and with weak Si and S lines) ; and type la 
SNe (SNe la), which show strong Si and S, but no H, lines in 
their spectra (see Filippenko ||1997 for a review; see Perets 
et al.|20To for a possible exception). CC SNe occur in mas- 



sive stars that have reached the end of their fuel cycles. 
Pre-explosion images have revealed directly the progenitors 
of some CC SNe, confirming that the progenitors of SNe II-P 
and SNe Iln are red and blue supergiants (or luminous blue 
variables), respectively; that most SNe Ib/c are the result 
of moderate-mass interacting binaries; and that broad-lined 
SNe Ic are the explosions of massive Wolf-Rayet stars (see 
ISmartt 12009] for a review). 

In contrast, SNe la are thought to be the result of the 
thermonuclear combustion of carbon-oxygen white dwarfs 
(WDs) that approach the Chandrasekhar limit through 
mass accretion in close binary systems (see Hillebrandt & 
|Niemeyer|[2000| and Howell | |2010] for reviews). Two basic 
routes have been suggested for the WD to grow in mass. 
The single-degenerate model postulates mass accretion from 
a main-sequence or giant companion star ( [Whelan fc Iben] 



1973), whereas the double-degenerate (DD) model invokes 



the merger of two WDs (Iben & Tutukov 1984 Webbink 
|1984[ ). However, there have been no unambiguous identifica- 
tions of SN la progenitors in pre-explosion images, or of re- 
maining companions in historical SN l a remnants (e.g.,|Voss 
^Nelcmans 2008; R oelofs et al.||2008| [Gonzalez Hernandez 
et al.|j2009i [Kerzendorf et al.||2009[ ). Programmes that seek 



to determine the DD merger rate by surveying for WD bina- 
ries (|Napiwotzki et al.|2004|[Geier et al.|2007[|Badenes et al.| 



2009 1 have yet to conclude whether this channel can account 



for some or all of the SN la rate. Thompson (20101 has re- 
cently proposed that at least some of the SN la progenitors 
may be triple systems, comprised of a WD- WD inner binary 



and a tertiary that induces Kozai ( 1962 1 oscillations in the 



inner binary, driving it to higher eccentricity and shortening 
the time until a gravitational-wave-driven merger between 
the two WDs. The possibility of detecting such triple sys- 
tems through their gravitational- wave signals is explored by 



Gould (20111 



One way to constrain indirectly the different SN la 
progenitor models is through their delay-time distribution 
(DTD) — the distribution of times between a hypothetical 
(5-function-like burst of star formation, and the subsequent 
SN la explosions. Different progenitor and explosion mod- 



els predict different forms of the DTD (e.g., 


Yungelson & 


Livio||2000| |Han & Podsiadlowski||2004| |Ruiter, Belczynski, 


& Fryer 


2009 


Mennekens et al. 


|2010|). Metallicity effects 



can also affect the DTD in some models (e.g., Kobayashi 



fc Nomoto 200"9] ). There are various ways to estimate the 



DTD observationally. Mannucci et al. ( 2005 \ compared the 



SN la rate per unit mass in different types of galaxies and 
found that the rate in blue galaxies is a factor of 30 larger 
than in red galaxies. This result led to the so-called 'A + B' 



model ( Scannapieco fc Bildsten|2005 I , which reproduces the 
SN la rate with a term proportional (through A) to the total 
stellar mass of the SN host population, and a second term 
which is proportional (through B) to the star formation rate 
(SFR) of the host population. The A + B model is effectively 
a two-time-bin approximation of the DTD. 

Totani et al. ( 2008[ ) compared the SN rates in ellip- 



tical galaxies in the Subaru-XMM Deep Field (SXDF) to 
the mean ages of their stellar populations, and deduced a 
power-law shape of the form ^'(t) oc t'^ for the DTD, with 
P ~ —1 in the delay-time range of 0.1-4 Gyr. [Maoz et al.| 



(20111 compared the SN rate and the star formation histo- 
ries (SFHs) of a subset of the galax;ies monitored by the Lick 
Observatory SN Search ( [Leaman et al.|[2011| |. They recon- 
structed a falling DTD, with significant detections of both 
'prompt' SNe la (with delays of < 420 Myr) and 'delayed' 
ones (> 2.4 Gyr). Similar results were obtained by [Brandt] 



et al. (20101, analysing the SNe la from the Sloan Digital 
Sky Su rvey II (SDSS-II; |York et al.|2000] ). |Maoz fc Badenes 



( [20101 ) compared between the SN rate in the Magellanic 
Clouds as inferred from SN remnants and the SFHs of their 
resolved stellar populations, and detected a prompt compo- 
nent in the DTD. Comparisons of the rates of SNe la and 
the luminosity- weighted mean ages of their host populations 



have been undertaken by Aubourg et al. ()2008 1; Raskin et al. 
(2009); Coope r, Newman, fc Yan| ( [2009[ ); [Sch awinski '(2009[): 
and Yasuda fc Fukugita (20101. While some of these stud- 



ies may be susceptible to biases resulting from the choices of 
'control samples' (see, e.g., Maoz [2008 \, they have generally 
also found evidence for a population of SNe la with short 
delays. 

Measurement of SN rates versus redshift in galaxy clus- 
ters has provided another powerful probe of the DTD. Clus- 
ter SFHs are relatively simple, and thus the form of the DTD 
is obtainable almost directly from the SN rate as a function 
of cosmic time. Furthermore, the deep gravitational poten- 
tials mean that the total metal content of clusters, as quan- 
tified by optical and X-ray measurements, provide a record 
of the time-integrated contributions, and hence numbers, of 
SNe over the cluster histories. This sets the integral of the 



DTD. [Maoz, Sharon, fc Gal- Yam (20101 have recently com- 



piled and analysed cluster SN rates from a number of sur- 
veys in the redshift range < z < 1.5 (Gal- Yam, Maoz, fc 



Sharon 2002; Sharon et al. 2007, 2010; Graham et al.||2008| 



Ma nnucci et al., 2008, Dilday et al. 2010b , Barbary et al.| 
' 2010| . They find that the best-fitting DTD is a power law 
with an index oi P = -1.1 ±0.2 or /3 = -1.3 ±0.2, depend- 
ing on the assumed value of the present-day stellar-to-iron 
mass ratio in clusters. Thus, a variety of recent attempts to 
recover the DTD, involving a range of techniques, redshifts, 
and environments, consistently indicate a power-law DTD 
with index /3 « — 1 (see Maoz et al.]|2010 for an intercom- 
parison of these results) . 

There is, however, one approach to recover the DTD 
that has produced some conflicting results. The SN rate in 
field galaxies at cosmic time t, i?ia(t), is the convolution of 
the SFH, S{t), with the DTD, *(t): 



i?ia(t) = / 5(t-r)*(r)dr. 
Jo 



(1) 



The DTD can therefore be recovered, in principle, by com- 
paring the field SN la rate vs. redshift to the cosmic SFH. 



Supernovae in the Subaru Deep Field 3 



The cosmic SFH has been measured out to z « 6 (see, 
e.g., the compilation of [Hopkins fc Beacom|[2b06[ hereafter 
HB06), and several surveys have attempted to extend these 
measurements out to z ~ 8 (|Verma et al. 120071 [Yiiksel et al. 



2008 



2009 



hereafter Y08;|Bouwens et al ||2 008 Rod dy fc Steidel 



Kistler et al.|2009| |Yan et al.|200 9). While all surveys 
observe a rise in the SFH towards z — 1-2.5, to date esti- 
mates of the SFH based on the ultraviolet (UV) emission 



of field galaxies (e.g., Bouwens et al.||2010J have produced 
shallower evolutions than those based on the far-infrared 
(IR) continuum of galaxies, (e.g., Le Floc'h et al.|2005 Ru- 



jopakarn et al.pOlO l. This is due to the systematic uncer- 
tainty introduced by the need to correct the observed UV 
luminosity for extinction by dust. A recent attempt by |Qda| 



et al. (2008 hereafter O08) to derive the cosmic SFH us- 
ing CC SN and SN la rate measurements found constraints 
which are consistent with the latest IR-based SFH measure- 
ments, and slightly higher than the latest UV-based mea- 
surements. 



Gal- Yam & Maoz ( 2004 1 were the first to set constraints 
on the DTD with this approach, based on a small sample of 
SNe la out to z = 0.8. A number of surveys over the past 
decade have measured the SN la rate out to 2 « 0.2 ( Cappef 



laro, Evans, & Turatto|1999| [Hardin et al. 


20001 |Pain et al. 


2002 


Tonry et al.||2003| Blanc et al.||2004 


[Botticella et al. 


2008 


hereafter B08;'Horesh ct al.'2008; Li ct al.|2011b|). Ad- 



ditional surveys, such as the SDSS (Madgwick et al.||2003 
[Dilday et al.|2008 |2010a I and the Supernova Legacy Survey 



(SNLS; |Neill et al.|2006| hereafter N06 |Neill et al.|2007| ) have 
added measurements out to z 0.8. The previously dis- 
cordant measurements of the Institute for Astronomy (IfA) 
Deep Survey ( Barris fc Tonry|2006 \ have recently been cor- 
rected and extended to redshift z = 1.05 by |Rodney fc Tonry| 



(20101 



Measurements of the SN rate at z > 1 were first real- 



ized by Dahlen et al. (2004 hereafter D04), using the Hub- 
ble Space Telescope {HSTj Advanced Camera for Surveys 
(ACS) observations of the GOODS fields. Additional data 
were analysed by Dahlen, Strolger, fc Riess (2008 here- 



after DOS). D04 and D08 argued that their data indicate 
a peak in the SN rate at z ~ 0.8, with a steep decline 
at higher redshifts. Based on this rate evolution, [Strolger] 



et al. (20041, D04, and DOS deduced a best-fitting narrow 



Gaussian-shaped DTD, centred at a delay time of 3.4 Gyr. 



Similarly, Strolger, Dahlen, fc Riess ( |2010) adopted a uni- 
modal, skew-normal function (see their equation 6) for the 
DTD, from which they inferred that the DTD should be con- 
fined to a delay-time range of 3-4 Gyr. However, analysing 
much of the same data, Kuznetsova et al. (|200S) found that 



they could not distinguish between a flat SN rate at z > 0.5 
and a decline at z > 1, due to the large statistical and sys- 
tematic uncertainties in the /fST/GOODS dataset. 



Horiuchi fc Beacom (20101 recently found that when 



coupled with the Y08 SFH, the Gaussian DTD proposed by 
DOS, along with the bimodal DTD from [Mannucci, Delia] 



Valle, fc Panagia (20061, underpredicted precise SN la rate 
measurements a.t z < 0.3. A power-law DTD with index 
j3 = — 1.0±0.3, however, fit the data well. A similar attempt 
by Blanc fc Greggio (2008) to couple between the cosmic 
SFH and the SN la rates from the above data also led to 
the conclusion that a broad range of DTD models could be 
accomodated by the data, including power-law DTDs, due 



to small-number statistics. In three HST cycles, GOODS 
found 53 SNe la, of which only 3 were in the lA < z < 
1.8 range. Larger SN la samples are clearly needed in order 
to determine precise rates at these redshifts, to recover the 
DTD, and to compare it to other measurements. 

To address this problem, in 2005 we initiated a ground- 
based high-z SN survey with the objective of determining 
the SN la rate at 2 > 1. Our survey is based on single-epoch 
discovery and classification of SNe in the Subaru Deep Field 



(SDF; [Kashikawa et aL][2004[ hereafter K04). In 2007 we 
presented initial results from our survey for SNe la out to 
z = 1.6, based on the first epoch of observations (Poznanski 



et al. 2007b hereafter P07b). This first epoch produced a 
number of SNe la that was similar to that found by D04 
in GOODS. The high-2 rates we found were also consistent 
with those of D04 and DOS, with similar uncertainties, but 
our results suggested a flat rather than a declining SN la 
rate at high redshifts. 

In this paper, we present our final sample of 150 SNe, 
based on four SDF epochs, and derive the most precise SN la 
rates to date at 1 < 2 < 2. In Section [2] we describe our ob- 
servations of the SDF and spectroscopy of several of our SN 
host galaxies. Sections [3] and [4] detail our methods for dis- 
covering the SNe and their host galaxies. In Section [5] we 
classify the SN candidates into SNe la and CC SNe with the 
SN Automatic Bayesian Classifier (SNABC) algorithm of 



Poznanski, Maoz, fc Gal- Yam (2007a) hereafter P07a). The 



distribution of SNe among types and redshift bins is exam- 
ined in Section [6] and corrected for biases introduced by the 
SNABC. We derive the SN la and CC SN rates in Section[7l 
The SN la rates, along with rates collected from the litera- 
ture, are then used to constrain the DTD in Section [S] The 
best-fitting DTD is used to predict the SN la rate at z > 2 
and calculate the accumulation of iron in the Universe, as a 
function of redshift, in Section|9] We summarise and discuss 
our results in Section [lO] Throughout this paper we assume 
a A-cold-dark-matter (ACDM) cosmological model with pa- 
rameters f^A = 0.7, ilm = 0.3, and Ho — 70 km s~^ Mpc~^. 
Unless noted otherwise, all magnitudes are on the AB sys- 
tem ( jOke fc Gun"n|T983| ). 



2 OBSERVATIONS AND REDUCTIONS 
2.1 Imaging 

The SDF {a = 13''24™39", 5 = -h27°29'26"; J2000) was first 
imaged by K04 with the Suprime-Cam camera on the Sub- 
aru 8.2-m telescope on Mauna Kea, Hawaii. Suprime-Cam 
is a 5 X 2 mosaic of 2k x 4k pixel CCDs at the prime focus 
of the telescope, with a field of vie w of 34 x 27 arcmin^ , 
and a scale of 0.202 arcsec pixel"^ (Miyazaki et al. 2002 1. 



K04 imaged the SDF in five broad-band filters {B, V, R, i , 
and 2') and two narrow-band filters {NB816 and NB921), 
over an area of 30 x 37 arcmin^, down to 3a limiting mag- 
nitudes oi B = 28.45, V = 27.74, R = 27.S0, i' = 27.43, 
2' = 26.62, NB816 = 26.63, and NB921 = 26.54 (5cr fimits 
of B = 27.87, y = 27.15, R = 27.24, i' = 27.01, 2' = 26.06, 
NB816 = 26.24, and NB921 = 26.07), as measured in circu- 
lar apertures having radii of 1 arcsec. See K04 for details of 
those observations. This initial epoch of optical observations 
is denoted here as 'epoch 1.' 
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Table 1. Summary of optical imaging data for epochs 2 through 5 



Epoch 


Band 


Exp. 


Seeing 


3(T miim " 


5(T miim * 


mo " 


UT Date 






¥\ 


[arcsec] 


[mag] 


[mag] 


[mag/count] 




2 


R 


7,920 


1.06 


27.18 


26.63 


33.93 


2005 Mar. 5/6 




i' 


10,800 


0.99 


27.00 


26.45 


33.99 


2005 Mar. 5/6 






18,240 


1.03 


26.33 


25.77 


32.92 


2005 Mar. 5/6 


3 


R 


11,460 


0.79 


27.98 


27.43 


34.08 


2007 Feb. 12/13/14/15 




i' 


15,000 


0.80 


27.79 


27.24 


34.11 


2007 Feb. 12/13/14/15 




z' 


27,240 


0.85 


26.90 


26.35 


33.01 


2007 Feb. 12/13/14/15 


4 


R 


8,220 


0.90 


27.36 


26.80 


33.14 


2007 May 15/16 




i' 


7,960 


0.84 


27.17 


26.62 


33.16 


2007 May 15/16 






17,150 


0.73 


26.86 


26.30 


31.87 


2007 May 15/16 


5 


R 


10,550 


0.83 


27.70 


27.14 


34.00 


2008 Jun. 1/2/3/4 




i' 


12,960 


0.81 


27.50 


26.94 


34.06 


2008 Jun. 1/2/3/4 




z! 


23,500 


0.73 


27.21 


26.66 


32.99 


2008 Jun. 1/2/3/4 



"3(7 hmiting magnitude, within a circular aperture having a radius the size of the image's seeing FWHM. 
^ba limiting magnitude. 

'^Magnitude zero point, i.e., the magnitude of a source in the image with 1 count (2.6 e^). 



In our analysis, we also make use of additional existing 
data on the SDF, particularly for estimating the proper- 
ties of the galaxies hosting the SNe we find. Near-infrared 
(NIR) photometry, in J and K, was obtained with the Wide- 
Field Camera on the United Kingdom Infrared Telescope 
(UKIRT; [Hayashi et al.||2009| Motohara et al., in prepa- 
ration) down to 3(T limiting magnitudes of J = 24.67 and 
K — 25.07 in apertures with radii of 1 arcsec {5a limits of 
J = 24.33 and K = 24.52 mag). While the /f-band data 
cover the same area of the SDF as the optical observations, 
the J-band data cover only ~ 40 per cent of the field. UV 
observations of the SDF were obtained by the Galaxy Evo- 
lution Explorer {GALEX; Ly et al.]|2009 1, with total expo- 
sures of 81 ks in the far-UV (FUV) band (A ^ 1530 A) and 
161 ks in the near-UV {NUV) band (A ^ 2270 A). These 
integration times result in 3cr limiting magnitudes of 26.42 
and 26.49 in the {FUV) and {NUV) bands, respectively, in 
apertures with radii of 7.5 arcsec (or 5a limits of 25.86 and 
25.93 mag). 

We reimaged the field on four separate epochs (UT 
dates are used throughout this paper): 2005 March 5 and 
6 (epoch 2, analysed by P07b); 2007 February 12-15 (epoch 
3); 2007 May 15 and 16 (epoch 4); and 2008 June 1-4 (epoch 
5). During epochs consisting of two nights, the SDF was ob- 
served during most of the night. On the epochs that were 
spread over four nights, either the first or the second half of 
each night was dedicated to the SDF programme. In either 
case, we consider the consecutive nights to be a single epoch, 
whose nightly images can be coadded, given the longer time 
scales on which SNe evolve at the redshifts we probe. On all 
occasions, we imaged the field in the three reddest Suprime- 
Cam broad bands: R, i' , and z' . These filters, which probe 
the rest-frame blue emission of SNe a,t z — 1-2, are the most 
suitable for discovering and classifying such SNe (e.g., |Poz-| 



nanski et aL]|2002[ |Gal-Yam et al.||2004| |Riess et al.]|2004 1 

We followed a dithering pattern similar to the one described 
by K04. Table [l] lists the exposure times, average seeing, 
and limiting magnitudes in each band, for epochs 2 through 
5. In general, the average seeing for each night ranged be- 
tween 0.7 and 1 arcsec full width at half-maximum intensity 
(FWHM). 



We reduced the Subaru observations with the Suprime- 
Cam pipeline SDFRED (Yagi et a l.|2002||Ouchi et al.|2004 |. 
Briefly, the individual frames were overscan subtracted, flat 
fielded using superflats, distortion corrected, sky subtracted, 
registered, and combined. In contrast to K04 and P07b, we 
did not apply point-spread function (PSF) degradation on 
the new images, since it reduces the frame depth. The com- 
bined image was then matched to the i'-band image from 
K04 by using the astrometri f] code to fin d the astrometric 
correction, and the IRAF|^ ( |Tody 1986 1 task wregister to 
register the two images. The photometric calibration of the 
images from epoch 1 was done by K04, achieving a precision 
for the zero points of ~ 0.05 mag (see section 4.2 of K04). We 
calibrated our images relative to epoch 1 by comparing the 
photometry of all the objects detected with SExtractor 
( Berlin &: Arnoutsl|1996 1 in both epochs. The mean of the 
differences between the two measurements was taken to be 
the difference in zero points. 

In order to create a reference image to be compared 
to each epoch, the images of all the other epochs were 
scaled, weighted according to their depth, and stacked using 
the IRAF task imcombine. The stacking process included 
a sigma-clipping procedure that excluded any transient or 
highly variable objects from the resulting summed image. 
Four 'master' images were created in this fashion, for each 
search epoch, where each such image is composed of all other 
epochs, except the search epoch in question. These master 
images proved deeper and sharper than the original epoch- 1 
images used by P07b as reference images for the subtraction 
process. For example, the epoch-5 master image has a So- 
limiting magnitude of ^ — 27.01, as measured in an aper- 
ture having a radius the size of the image's PSF FWHM 
of 0.96 arcsec, and is the deepest of the master images. As 
discussed in Section |3.1| below, the use of the new master 



^ http:/ /www. na.astro.it/~radovich 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation (NSF). 
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images as reference images resulted in the discovery of SNe 
in epoch 2 that went undiscovered by P07b. 

We performed PSF matching, scaling, and image sub- 
traction between the target and reference images in each 
Subaru epoch in all bands, using the software HOTPANT^ 
an implementation of the ISIS algorithm of |Alard fc Lup 



The DEIMOS data were reduced usiiig a modified ver- 
sion of the DEEP2 data-reduction pipelina*! which bias cor- 
rects, flattens, rectifies, and sky subtracts the data before 
extracting a spectrum (Foley et al. 20071. The wavelength 



Becker 



ton| (|1998|l for image subtraction (as described by 
et al . 2004). Briefly, HOTPANTS divides the images into a 
predetermined number of regions, and in each region finds 
the convolution kernel necessary to match the PSF of one 
image to that of the other. HOTPANTS is similar to ISIS, 
which was used by P07b, but allows more control over the 
subtraction process. For example, each region of the image 
is subdivided into stamps and substamps, where the sub- 
stamps are centred on astronomical objects. The kernel is 
then computed for each substamp, producing a distribution 
of values used to sigma-clip outliers, thus ensuring a more 
accurate determination of the kernel in each stamp, and ul- 
timately a better mapping of the spatial variations of the 
kernel across the image. We also made use of the software's 
ability to mask saturated pixels, which vastly reduced the 
number of residuals in the difference images. 

As a consequence of the dithering, the final images have 
a field of view of 0.31 deg'^; however, due to the different ef- 
fective exposures in the fringes of the field, a substantial re- 
gion along the edges suffers from a significantly lower signal- 
to-noise ratio (S/N). We therefore crop the edges of the dif- 
ference image, ending with a total subtraction area of 0.25 
deg^ 



2.2 Spectroscopy 

As detailed in section 2.2 of P07b, we obtained spectra of 
17 of the SN host galaxies from epoch 2, together with sev- 
eral hundred random galaxies in the SDF, using the Low- 



Resolution Imaging Spectrometer (LRIS; Oke et al. 1995 1 
on the Keck I 10-m telescope, and the Deep Imaging Multi- 
Object Spectrograph (DEIMOS; |Faber et alT||2003| on the 
Keck II 10-m telescope. 

In addition to the SN host spectra published by P07b, 
we obtained spectra of 7 additional SN host galaxies. These 
spectra were taken during observations carried out on the 
night of 2010 February 15 with DEIMOS on the Keck II 
telescope. The single mask utilised for these observations 
contained 16 SN host galaxies, as well as the positions of 
tens of filler galaxies. The mask was observed for a total 
of 3 X 30 min. We used the 600 line mm~^ grating, with 
the GG455 order-blocking filter and a wavelength range of 
~ 4400-9600 A, with the exact limits depending on each 
individual spectrum. 

The 600 line mm"'^ grating yields a FWHM intensity 
resolution of ~ 3 A, or ~ 120 km s~^, at 7500 A. This res- 
olution is sufficient to resolve many night-sky lines and the 
[O II] AA3726, 3729 doublet. By resolving night-sky lines, 
one can find emission lines in the reddest part of the spec- 
trum, where sky lines are blended in low-resolution spectra. 
Furthermore, by resolving the [O II] doublet, we can confi- 
dently identify an object's redshift, even with only a single 
line. 



solutions were derived by low-order polynomial fits to the 
lamp spectral lines, and shifted to match night-sky lines at 
the positions of the objects. Finally, the spectra were fiux 
calibrated by scaling them to the mean fiuxes in the R and 
/ bands. Consequently, the displayed continuum spectral 
shape is not precisely calibrated. In any event, the contin- 
uum emission of the host galaxies is weak and noisy, and 
therefore we rely on spectral lines alone for redshift identi- 
fication. 



3 SUPERNOVA CANDIDATES 

In this section we describe the methods by which we have 
discovered the SN candidates in our sample, derive the de- 
tection efficiency of the survey, and measure the photomet- 
ric and astrometric properties of the candidates and their 
uncertainties. We have discovered a total of 163 transient 
objects, of which 150 are most likely SNe. The luminosities 
of the transients, inferred from their measured photometry 
and the redshifts of their associated host galaxies (as de- 
below). 



rived in Section 4.2 



lead us to conclude that these 
150 events are SNe. In Section 13.11 we describe the crite- 
ria according to which the transients were chosen, culling 
random noise peaks, image subtraction artefacts, and previ- 
ously known active galactic nuclei (AGNs). We calculate the 
probability of contamination by flaring Galactic M dwarfs 
and unknown AGNs in Section |4.1| The probable contam- 
ination by AGNs is compared with the number of actual 
possible AGNs among the candidates in Section [5. 1| In Sec- 
tion |4.1| we also calculate the probability of a chance associa- 
tion between a transient object and its surrounding galaxies. 

Since our survey classifles SNe based on single-epoch ob- 
servations without spectroscopic follow-up observations, the 
SNe we discover do not satisfy the International Astronomi- 
cal Union's criteria for a 'standard' designation. As in P07b, 
we will continue to use the following naming conventions. 
We denote the SNe from epochs 2 through 5 respectively 
as 'SNSDF0503.XX,' 'SNSDF0702.XX,' 'SNSDF0705.XX,' 
and 'SNSDF0806.XX,' with the first two digits denoting 
the year, the next two digits the month, and XX being a 
serial number ordered according to the SN z'-band appar- 
ent magnitude. The respective host galaxies are referred to 
as 'hSDF0503.XX,' 'hSDF0702.XX,' 'hSDF0705.XX,' and 
'hSDF0806.XX.' 



3.1 Candidate selection 

The a'-band difference image obtained with HOTPANTS 
was scanned with SExtractor to search for variable ob- 
jects. SExTRACTOR was set to identify and extract all objects 
which had at least 6 connected pixels with flux 3a above the 
local background level. T. Morokuma (private communica- 
tion) provided us with a catalogue of 481 AGNs, which were 
identified in epoch 1 by their long-term i'-band variability. In 



^ http:/ /www. astro. washington.edu/users/becker/hotpants. html 



^ http: / / astro. berkeley.edu/'~cooper/deep/spec2d/ 
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Figure 1. SN host-galaxy spectra from the 2010 February 15 Keck DEIMOS observations, with the prominent emission and absorption 
features marked. The spectra have been rebinned into 10 A bins, (a) hSDF0702.03, z = 0.70; (b) hSDF0702.21, z = 0.30; (c) hSDF0702.23, 
z = 0.96; (d) hSDF0705.18, z = 1.41; (e) hSDF0806.48, z = 1.13; (f) hSDF0806.54, z = 0.53; and (g) hSDF0806.55, z = 0.60. 



our survey, these galaxies were therefore ignored, as further 
discussed in Section [4.1| These galaxies still constitute fewer 
than 1 per cent of all galaxies in the SDF, and therefore this 
exclusion has negligible effect on our SN survey. 

In order to reject other non-SN events, the remaining 
variable candidates were examined as follows. 

(i) Of the objects identified by SExtractor, we rejected 
all those which showed suspect residual shapes, indicative 
of a subtraction artefact. For maximum completeness, the 
threshold for SExtractor detection was set low, and thou- 
sands of candidates were inspected by eye by one of us (OG). 

(ii) We compared two 2;'-band difference images of the same 
field. The main difference image was obtained by allowing 
HOTPANTS to calculate the convolution kernel for the sub- 



traction over the entire image. A second difference image 
was obtained by forcing HOTPANTS to break the image 
into four subregions, and calculate the convolution kernel 
in each one. This second difference image was generally less 
clean than the first, but allowed for the rejection of subtrac- 
tion artefacts in the main difference image, as not all of those 
would be reproduced in the second subtraction process. 

(iii) We compared the main 2;'-band difference image in a 
certain epoch with difference images of the other epochs in 
order to identify and reject AGNs that were not already 
rejected based on the Morokuma AGN catalogue, or other 
objects that exhibited variability over a large stretch of time. 
Roughly 40 transients were identified as AGN candidates 
due to their variability over several epochs. These objects 
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were not included in the Morokuma AGN catalogue, and 
may have been quiescent at the time it was compiled. 

(iv) In order to further reject subtraction artefacts, we 
stacked the exposures in each epoch into two subepoch im- 
ages, where each subepoch was composed of half of the ob- 
servation nights. These images were then used to obtain new 
difference images which we compared with the main ^'-band 
images. As in the previous steps, objects which appeared in 
the main difference image, but not in the subepoch differ- 
ence images, were rejected. We note that solar-system ob- 
jects were already eliminated in the nightly averaging, since 
even as far as 30 AU (Stern & Colwcll 1997. ) a Kuiper Belt 
object would move due to the Earth's motion by ~ 40 arcsec, 
or 200 pixels, in the course of an 8-hour night. 

(v) For every candidate found in the 2^ band, difference im- 
ages in the R and i' bands were also examined, and objects 
which showed suspect residual shapes were rejected. We note 
that no candidate was rejected because of a nondetection in 
the R or i' bands, since at least some high-z SNe are ex- 
pected to be very faint or undetected in the observed-frame 
R and i' bands. 

(vi) Finally, for each SN candidate, we derived the local S/N 
by dividing the SN counts in an aperture of I arcsec radius 
(before application of an aperture correction) by the stan- 
dard deviation of the total counts in tens of identical aper- 
tures centred on surrounding blank regions. SN candidates 
which had a S/N smaller than 3 were rejected as probable 
noise peaks. 

We note that steps (ii) and (iii) are selection criteria addi- 
tional to those followed by P07b. 

In order to apply our new criteria uniformly to the full 
SN survey, we resurveyed epoch 2. Of the 33 SNe found by 
P07b, 28 were recovered. The SN candidates listed in P07b 
as SNSDF0503.27, SNSDF0503.33, and SNSDF0503.40 were 
not detected by SExtractor, because the S/N was too 
low. While the first two SN candidates listed above ap- 
pear in the difference images, we do not detect the third 
one in our renewed analysis. SNSDF0503.29 was detected 
by SExtractor, but whereas in the main difference im- 
age it appears as a point source, in the secondary differ- 
ence image it is extended, and the position of its centre is 
offset by ~ 0.35 arcsec. SNSDF0503.32 was not detected 
by SExtractor, and while it appears in the main differ- 
ence image, it is absent from the secondary difference im- 
age. Thus, with our improved reference images and image 
subtraction procedures, these events from P07b do not pass 
our current selection criteria. 

On the other hand, we have discovered 8 new SN 
candidates in epoch 2, not reported by P07b. In this 
work, these SN candidates are hsted as SNSDF0503.06, 
SNSDF0503.I6, SNSDF0503.I9, SNSDF0503.27, 
SNSDF0503.3I, SNSDF0503.32, SNSDF0503.33, and 
SNSDF0503.34. The differences between the P07b sample 
and the present sample are due to two reasons: (a) the 
use of HOTPANTS in the current work, which provides 
cleaner subtractions than ISIS, and (b) the use of deeper 
z'-band master images with better image quality, instead of 
the shallower epoch-1 ^/-band image, as references. In any 
event, the list of epoch-2 SNe that we report in Table [8] 
supersedes the one presented by P07b. 



3.2 Detection efficiency simulation 

In our survey, SNe may be missed as a result of many effects, 
including imperfect subtractions, noise fluctuations, and hu- 
man error. In order to quantify these systematic effects, we 
measure our detection efficiency by blindly planting artificial 
point sources, which match the SN population in our survey 
as closely as possible, in the presubtraction ^/-band images, 
and then discovering them along with the real SNe. The sim- 
ulated SN sample was constructed as detailed in section 3.2 
of P07b. Our resulting efficiency as a function of magnitude, 
in each epoch, can be seen in Fig. [2] We follow [Sharon et al.| 
( 2007 1 and fit the following function to the data: 



rj{m-Tni/2,si,S2) = 



'1/2 



M/2 



m ^ mi/2 
m > mi/2. 



(2) 

where m is the ^/-band magnitude of the fake SNe, mi/2 is 
the magnitude at which the efficiency drops to 0.5, and si 
and S2 determine the range over which the efficiency drops 
from I to 0.5, and from 0.5 to 0, respectively. 



3.3 Supernova sample 

We have found a total of 150 SNe, with magnitudes in the 
range z = 22.9 to z' = 26.7. Table [s] lists the SNe and 
their properties. Apart from these 150 SNe, we detect sev- 
eral tens of candidates at fainter magnitudes, as we would 
expect based on our efficiency simulations, but these are all 
objects with S/N < 3. While some of these objects may 
be SNe, an unknown number of them could be false posi- 
tives, such as subtraction artefacts or random noise peaks. 
We therefore limit our sample to 2/ < 26.6, zf < 26.4, and 
z' < 26.7 mag for epochs 3 through 5, respectively. These 
are the values of mi/2 in each epoch. In epoch 2 we reach 
the 50 per cent efficiency mark at 26.2 mag. However, in the 
interest of backward compatibility with P07b, we lowered 
the efficiency cutoff for epoch 2 to 26.3 mag. 

Using SExtractor, we have performed aperture pho- 
tometry of the SNe in the R, i' and 2/ difference images 
within fixed f-arcsec-radius circular apertures. To estimate 
the aperture correction and photometric uncertainty, we 
measured the magnitudes of ~ 600 simulated point sources, 
ranging in brightness from 23 to 28 mag, planted in a 4k x 
4k pixel subframe of the SDF R-, i'-, and z'-band images. 
We took the difference between the average of the magni- 
tude in each bin and the true magnitude as the required 
aperture correction, and the standard deviation in each mag- 
nitude bin to be the minimum photometric statistical error 
for objects of that magnitude. For example, the mean aper- 
ture correction for the epoch-2 z'-band image was 0.2 mag 
(i.e., due to aperture losses, the measured photometry was 
0.2 mag too faint) and the standard deviation ranged from 
0.03 to 0.29 mag from the brightest to the faintest artificial 
sources, respectively. The adopted uncertainty for each SN 
was taken to be the larger among the uncertainty computed 
by SExtractor and the statistical uncertainty for the given 
magnitude bin from the simulations. 

We also measured the offset of each SN from its host 
galaxy. To estimate the uncertainty of the offset, ~ 12,000 
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Figure 2. Fraction of simulated SNe recovered as a function of z'-band magnitude. Error bars indicate la binomial uncertainties. The 
dotted lines mark the 50 per cent efficiency mark. 



simulated point sources, divided into magnitude bins of 
width 0.3 mag, were planted in the z'-band image of 
each epoch. We then measured their locations, in both 
the original image and the ^'-band difference image, using 
SExTRACTOR, and took the mean of the location residuals 
in each bin as an estimate of the uncertainty of the object's 
location. This uncertainty was added, in quadrature, to the 
uncertainty in the location of the SN host galaxy. The real 
SN offsets ranged between and 3.61 arcsec, and the un- 
certainties ranged between 0.02 and 0.16 arcsec, with the 
centres of brighter sources being, of course, better localized. 



4 SUPERNOVA HOST GALAXIES 

In this section, we determine the host galaxy of each SN and 
then measure its properties. The SN host galaxies, including 
their photometry in all available bands, are presented in 
Table H 



4.1 Identification and pliotometry 

The SN host galaxies were chosen to be the closest galax- 
ies, in units of those galaxies' half-light radii, as measured 
with SExTRACTOR in the i' band. A small number of SNe 
had several possible hosts. To choose between them we mea- 
sured the photometric redshift (photo-z) of each host. If the 
different hosts were found to be at the same redshift, that 
redshift was adopted for the SN as well. If, on the other 



hand, the diflerent hosts were found to lie at different red- 
shifts, we computed the likelihood of a SN of the type, as 
classified by SNABC, at those different redshifts being ob- 
served at the magnitude measured. In this manner we were 
able to eliminate unlikely hosts. 



Using SExTRACTOR, we measured the Petrosian ( 1976 1 



magnitude of the host galaxies in the seven optical bands of 
epoch 1. We chose Petrosian photometry, since it measures 
the flux of resolved objects within a given fraction of the 
object's light profile, thus enabling one to compare between 
measurements taken in different filters. The resulting cata- 
logue was cross-matched with the J and K catalogues. Addi- 
tionally, for each galaxy we checked the corresponding loca- 
tion in the GALEX FUV and NUV background-subtracted 
images. Since the GALEX PSF is much larger than that of 
Subaru and UKIRT, most of our galaxies appear as point 
sources, making it impossible to measure Petrosian magni- 
tudes; hence, any measurement within any aperture would 
not capture the same percentage of light as in the optical and 
NIR bands. Furthermore, owing to the density of sources in 
the SDF and the size of the GALEX PSF, in many cases 
it proved impossible to determine which source in the opti- 
cal image was associated with the UV signal. In those cases 
where we could associate nondetections in the UV bands 
unambiguously with our host galaxies, we added the limit- 
ing magnitudes in the relevant UV bands to the catalogue. 
In Section |4.2| we detail how we combined these limiting 
magnitudes with the optical and NIR data to compute the 
redshifts of the SN host galaxies. As with the SN photom- 
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etry, for the host photometry we estimated the uncertainty 
in each magnitude bin using artificial sources with galac- 
tic profiles (created with the IRAF routine gallist) that we 
planted in the images. 

To test whether any of our chosen host galax;ies are 
merely chance associations, we counted the fractions of the 
total imaged SDF area that are within 0.1-light-radius-wide 
annuli of all the galaxies detected in the field. From this we 
conclude that, among the 110 SNe within 5C 0.5 light radii 
of their chosen hosts, < 1 SN is expected to be a chance 
association. These 110 SNe include all 12 SNe in the 1.5 < 
z < 2.0 range, and 24 of the 26 SNe in the 1.0 < 2 < 1.5 
range. At larger host-SN separations, 23, 6, and 1 of our 
SNe are found within 0.5-1.0, 1.0-1.5, and 1.5-2.0 hght radii 
of their host galaxies, respectively. Among these events, we 
expect 6, 3, and < 1 (respectively) to be chance associations. 
However, 28 of these 30 large-separation events are at z < 1. 
Thus, while some fraction of our 2 < 1 rate may be due to 
contamination by chance associations, we estimate that our 
1.0 < 2 < 1.5 rate is affected by such contamination at only 
the few-percent level, and the 1.5 < 2 < 2.0 negligibly so. 
In P07b we found that, assuming a |Sersic | ( [T968| 



model for the galaxy radial profile between n = 4, the |de| 
Vaucouleurs 1948 law (Peng et al. 20021 and n = 1, an 



exponential disk ( [Freeman |1970[|Peng et al.|2002[ ), between 
91 and 99.99 per cent of the light (respectively) falls within 
6 half-light radii of the galaxy's centre. Ten of our SNe 
have no visible host galaxies within this distance, and 
so we label them 'hostless' (namely SNSDF0503.14, 
SNSDF0503.18, SNSDF0702.06, SNSDF0705.20, 
SNSDF0705.21, SNSDF0705.24, SNSDF0806.04, 
SNSDF0806.30, SNSDF0806.49, and SNSDF0806.53). 
The probable host galaxy of SNSDF0806.51 appears 
exclusively in the B and R bands of epoch 1. Given that 
our photometric redshift estimate requires at least three 
photometric bands for its calculation, and that even the B 
and R detections are barely above the limiting magnitudes 
in those bands, we treat this SN as hostless as well. The 
most probable explanation is that these SNe occurred in 
galaxies fainter than the limiting magnitudes in all the 
photometric bands of epoch 1. 

Other possibilities to consider are that these candidates 
are high-2 AGNs or flaring Galactic M dwarfs. The fact that 
these hostless SN candidates are detected in only a single 
epoch over a period of 3 years argues against the AGN op- 
tion, as follows. Among the 481 objects identified in the 
Morokuma AGN catalogue, fewer than 1 per cent display 
detectable variability in only one of our four search-epoch 
difference images. 50 of the SN candidates in our sample lie 
within 0.2 arcsec (or 1 pixel) of their respective host-galaxy 
nuclei, and so could potentially be AGNs. Together with the 
above 11 hostless SNe, the predicted number of contaminat- 
ing AGNs in our sample is 61 x 0.01 « 0.6. The Poisson prob- 
ability of having at least one AGN in the sample is then ~ 45 
per cent, which is consistent with our having found one such 
object. The probability of finding two or more such objects 
drops to ~ 12 per cent (see S NSDF0 705.17 in Section [sTI 
and SNSDF0705.30 in Section IsTel). 



As to the second possibility, M-dwarf optical fiares con- 
sist of a fast rise followed by a decay lasting typically of 
order an hour or less, with the distribution of flare dura- 



2011 1. The longest known flares last ~ 10 hrs (Kowalski 
et al.|[2010 |, and these constitute < 1 per cent of all flares 
(E. Hilton, S. Hawley, private communication). With such 
variation timescales, M-star flare events would be flltered 
out in our nightly image averaging, or would at least show a 
decline between consecutive half-night averages. None of the 
hostless candidates show such a decline. We note, further, 
that flaring activity is limited to the younger M dwarfs in 
the Milky Way disk that are within a height of Z < 300 pc 
above the disk. Activity in older dwarfs, which have had time 



to be scattered to larger heights, is exceedingly rare (West 
etaLj2008; Ko walski et al . 2009; Walkowicz et al.|2011[ ). Any 



M dwarfs below the SDF detection limits in quiescence, and 
that had flared into visibility during our observations, would 
necessarily be at distances ^ 50 kpc, i.e., they would belong 
to the Galactic halo, and hence would be even older and less 
active than the Z > 300 pc disk stars. We therefore deem it 
highly unlikely that any of our hostless SN candidates are 
optical flares of Galactic M dwarfs. 

4.2 Host redshifts 

From our spectroscopy, detailed in Section |2.2[ we derived 
spectroscopic redshifts (spec-2) for 24 of the SN host galax- 
ies. Of these 24 SN host galaxies, hSDF0705.18 has the high- 
est spec-2, a.t z — 1.412. The seven new spectra obtained on 
2010 February 15 appear in Fig. [l] For the majority of our 
SN host galaxies, which are too faint for spectroscopy, we 
derive photometric redshifts, as in P07b, using the Zurich 



Extragalactic Bayesian Redshift Analyzer (ZEBRA; Feld- 
mann et al.|[2006l ). We calibrated ZEBRA in the manner 
described by P07b, but with a larger training set of 431 
galaxies, of which 150 are in the range 1 < 2 < 2. This train- 
ing set consisted of 123 galaxies imaged in the Keck runs 
detailed by P07b, along with data from other surveys that 
had been conducted in the SDF (e.g., Kashikawa et al.|2003 



tions steeply falling at longer durations (Walkowicz et al 



[20061 [Shimasaku et al.|2006l and a new sample obtained by 
N. Kashikawa in 2008 with DEIMOS on the Keck H tele- 
scope). ZEBRA was allowed to run over the redshift range 
< 2 < 3. 

Since ZEBRA does not, at the moment, offer an ad- 
equate treatment of upper limits, but rather deals with 
them as with any other photometry measurement, we de- 
cided (at the suggestion of R. Feldmann, private commu- 
nication) to halve the la FUV and NUV flux limits, and 
treat them as measurements with relative uncertainties of 
100 per cent, thus requiring ZEBRA's flt to pass through 
the region [0, /lo-]- If no GALEX signal existed that could 
be clearly associated with the optical galaxy, we used the 
UV flux limit (as described above) as an extra band in the 
ZEBRA flt, thus constraining the SEDs to those with fluxes 
lower than the UV flux limit. These upper limits on the UV 
flux were particularly useful for constraining the redshifts of 
galaxies having 'Lyman breaks' due to absorption by neu- 
tral hydrogen in the intergalactic medium (IGM). If, on the 
other hand, there was a GALEX detection, but due to the 
large GALEX PSF we could not clearly associate the UV 
signal with the optical SN host galaxy, we did not use the 
GALEX data at all. For larger samples, where more galaxies 
have clear signals in the UV, one could treat the UV signal as 
a lower limit, in similar fashion to our use of nondetections 
as upper limits. 
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Table 2. SN luminosity functions, presented as B-band absolute 
magnitudes (Vega) at maximum light, and Gaussian width. 
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Figure 3. Comparison of the photometric redshifts derived with 
ZEBRA and the corresponding spectroscopic redshifts for the 431 
galaxies in our training set (grey crosses) and for 24 SN host 
galaxies (empty diamonds). Error bars are the la confidence lim- 
its from the z-PDF of each galaxy. The rms scatter of the data is 
<7Az/(i+zs) = 0.065 for the training set and itaz/{i+Z(,) = 0.028 
for the SN host galaxies. 



Fig. |3] displays the ZEBRA photo-z values for our 
training-set galaxies, compared to their spec-z values. The 
training set of 431 galaxies has a root-mean square (rms) 
scatter of crAz/(l + Zs) = 0.075 (where Az — Zs — Zp) in the 
range < Zp < 2, after rejecting six 4(t outliers. This is con- 
sistent with the accuracy achieved by P07b, (JAz/(1 + Zs) — 
0.08, for 296 galaxies in the range < Zp < 1.8 and after 
rejection of five 4a outliers. The rms scatter for our 24 SN 
host galaxies is smaller: (taz/(1 + Zs) = 0.044. There were 
no 4(T outliers among these host galaxies. 

Of the various end products computed by ZEBRA, we 
use the redshift probability distribution function (z-PDF) of 
each SN host galaxy that results from marginalizing the full 
posterior distribution over all templates. In this manner the 
uncertainties in the determination of the photo-z are propa- 
gated into the classification stage. While most of the z-PDFs 
display a single, narrow peak, some are more structured, a 
result of degeneracies between the different combinations of 
redshifts and normalization constants (i.e., a certain galaxy 
may fit the same template if it is bright and distant, or if it 
is faint and nearby) or of a dearth of information. For exam- 
ple, the optical continuum shape of late-type galaxies can be 
approximated with a power law, and so its shape is weakly 
affected by redshift (see, e.g.. Fig. [7|. In such cases the UV 
data can be useful; a clear signal (whether a detection or 
a nondetection) in the NUV band would decide among the 
redshift values. In order to take the uncertainty introduced 
by the shape of the z-PDFs into account, we use the full 
2-PDFs in the classification stage (see Section [5|. 

For 23 of the 24 SN host galaxies with spectral red- 
shifts, the spec-z and photo-z values are almost identical, 
with A2/(l + z) < 0.08, while for hSDF0503.24 the differ- 
ence is only Az/(1 + z) = 0.10. For these galaxies we do 
not take the z-PDF computed by ZEBRA as input for the 
SNABC, but rather use a Gaussian z-PDF centred on that 
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galaxy's spec-z, with a width Wz = 0.01. For the 11 hostless 
SNe, we use a z-PDF which is the sum of the z-PDFs of all 
the other host galaxies. A different composite z-PDF, the 
average of the z-PDFs of all the galaxies in the SDF, was 
also tested for these SNe, and produced the same results. 
Given the resulting redshifts, the host galaxies of the host- 
less SNe would have to be fainter than between —15.8 and 
— 17.0 (absolute observed «'-band magnitude) to be unde- 
tected in the i'-band master images. This is consistent with 
these SNe having occured in low luminosity dwarf galaxies 
(see, e.g., Arcavi et al.||20l6 |. 



5 SUPERNOVA CLASSIFICATION 

We classify our SNe into SNe la and CC SNe using the 
SNABC algorithm of P07a. Briefly, the SNABC receives as 
input the photometry and z-PDF of a SN candidate. Using 
the above inputs, the SNABC then compares the colours of 
the SN candidate to the synthetic colours derived from a set 
of SN spectral templates of different types, ages, redshifts, 
host-galaxy and Galactic extinctions (based on the spectral 
templates of [Nugent, Kim, fc Perlmutter| |2002| hereafter 
N02)[^and to the rest-frame B-band luminosity functions 
(LFs) of the different SN types. In this work we used the LFs 
quoted by B08 for type la and II-P SNe, the LF measured 
by 



by 



Drout et al. (20101 for Ib/c SNe, and the LF measured 



Kiewe et al.| ( |2010| ) for type Iln SNe. [Drout et al.| ( |2010 l 

17.9 ± 0.9 for SNe lb 



measured peak magnitudes of Mr = 
and Mr = -18.3 ± 0.6 for SNe Ic. We take the weighted 
average of these magnitudes and get AIr — —18.2 ±0.9 
mag. Based on the N02 spectral template for SNe Ib/c, we 
apply a colour correction of {B ~ R) = 0.6 and arrive at 
A'Ib = —17.6 ± 0.9 for SNe Ib/c. In a similar vein, we apply 
a colour correction of {B — V) — —0.15 to the LF measured 
by Kiewe et al. ( |2010[ ) for SNe Iln, and arrive at a peak mag- 
nitude of Mb ~ —18.55 ± 1.00. The LFs and their sources 
are listed in Table[2] The host-galaxy extinction was allowed 
to vary in the range Av = 0-3 mag, which spans the full 
range of possible extinctions that we consider (see Section [6] 
for a discussion of the extinction model we use). 

The SNABC, as described by P07a, uses only the 
SN la and SN II-P spectral templates for classification. P07a 
describe how using more templates, such as SN Iln and 
SN Ib/c, allows for better classification of CC SNe, but at 
the same time significantly increases the number of SNe la 
misclassified as CC SNe, thus lowering the overall classifica- 
tion accuracy. We note that the goal of the current survey 
is not to discover and classify all types of SNe in the SDF, 



http: / / supernova.lbl.gov/~nugent / nugent_templates.html 
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Figure 4. ZEBRA fits and resultant redshift PDFs of tlie 1.5 < z < 2.0 SN la liost galaxies. The left panel of every pair shows the actual 
photometry (filled circles), the best-fitting galaxy template (solid line), and its synthetic photometry (empty circles). The vertical error 
bars denote the photometric uncertainty, and the horizontal error bars show the width of the filter. The header gives the designation of 
the SN host galaxy, most probable photo-z (zp), the spcc-z (zs, if such a measure exists for the specific object), the P^r degree of 
freedom of the fit, and the absolute B-band magnitude the object would have at Zp. The right panel of every pair shows the resultant 
z-PDF. If a spec-z exists for the SN host galaxy, it appears as a cross. 



but rather to determine the rates of SNe la statistically. 

The SNABC was designed and discussed specifically with 
the SDF survey, and its statistical approach, in mind. 

The SNABC computes the likelihood of each compari- 
son, and then marginalises over age, redshift, and extinction 
to arrive at the 'evidence' that the candidate is of a certain 



type: E{Ia) and E{CC). The evidence is then used to de- 
rive the probability that the candidate is either a SN la or 

CC SN, according to 

^^^^^ = £(Ia)Ts(CC)- 
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Figure 5. SNe la and host galaxies at 1.5 < z < 2.0. North is up and cast is left. All tiles are 10 arcsec on a side. The left-hand tiles 
show the SN host galaxies as imaged in epoch 1, whereas the centre tiles display the SN host galaxy as imaged in epochs 2 through 5. R-, 
i'-, and jZ-band images were combined to form the blue, green, and red channels (respectively) of the color composites. The right-hand 
tiles show the subtraction in the 2' band. Whereas the stretch of the colour images differs from panel to panel in order to highlight host 
properties, the greyscale for all difference images is identical. The header of each panel gives the designation of the SN la and its rcdshift. 
Similar images of the full sample of SNe are available in the electronic version of the paper. 



In addition to P(Ia), for each SN type the SNABC also 
produces a posterior z-PDF, which is constrained by the 
prior 2-PDF input from ZEBRA. The SNABC also produces 
a value that indicates how well the SN's colours compared 
with those of the best-fitting spectral template. A high 
value may imply the SN is a peculiar type of SN, an AGN, 
or a subtraction residual which was not rejected earlier. An 
event is considered a SN la if -P(Ia) > 0.5 (and a CC SN 
if -P(Ia) < 0.5). P07a have shown that P(Ia) can also be 
viewed as a confidence estimator: the closer it is to unity 
(zero) , the safer the classification of the candidate as a SN la 
(CC SN). P07a also found that for the sake of classification, 
most CC SNe resemble SNe II-P more than SNe la. Thus, 
while SN la classifications usually result in small values 
(X^ < 1), CC SN classifications may result in higher values, 
since SNe Iln or SNe Ib/c are forcibly compared to SN II-P 
spectral templates. 

The posterior redshift assigned to each SN by the 
SNABC usually matches the prior redshift assigned by ZE- 
BRA to within 5 per cent. In those cases where the difference 
between the two exceeds 5 per cent, we check the shape of 
the 2-PDF. A wide or multi-peaked z-PDF implies that the 
colours of the SN provided either more information than the 
a-PDF itself, or enough information to break the degeneracy 
between the different peaks in the z-PDF. In such instances 
(20 of the 150 SNe in our sample), we use the posterior red- 



shift computed by the SNABC. For example, SNSDF0806.32 
has a posterior redshift of 1.66, even though this value cor- 
responds to the weaker of the two peaks in the z-PDF of 
hSDF0806.32, as shown in Fig.[4] 

Table [8] lists the SNe in our sample, along with their 
redshifts and classifications. Of the 150 SNe in our sample, 
26 were found in the z < 0.5 bin, of which 5 were classified as 
SNe la and 21 as CC SNe. The 0.5 < 2 < 1.0 bin contains 86 
SNe, of which 50 were classified as SNe la and 36 as CC SNe. 
The 1.0 < z < 1.5 and 1.5 < z < 2.0 bins contain 26 and 
12 SNe, respectively, all of which were classified as SNe la. 
Two of the 12 SNe in the 1.5 < z < 2.0 bin have high 
values, and are dealt with individually in Section [5.1.6| The 
remaining 10 high-z SNe la are shown in Fig. [5] 

5.1 Notes on individual supernovae 

The high x^ values (> 10) of some of the 163 transients in 
our sample prompted their reevaluation and, in some cases, 
rejection. The final sample, after such rejections, includes 
150 SNe. All Xr values quoted are per degree of freedom. 

5.1.1 SNSDF0503.25 

While SNSDF0503.25 was classified as a CC SN [P(Ia) = 
0.06] with a high Xr value of 32, it is displaced from the 
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Figure 6. The possible host galaxies of SNSDF0702.30. The 
arrows in the x'-band image point to the two possible hosts: 
hSDF0702.30a is the resolved galaxy, while hSDF0702.30b is a 
compact source to the NW (above and to the right) of the latter. 
While there may be an ambiguous detection in the NUV band, 
both galaxies are clearly undetected in the FUV band. All tiles 
are 10 arcsec on a side. 



nucleus of its spiral host by 0.63 ± 0.07 arcsec. This, to- 
gether with the absence of the object at other epochs, ar- 
gues against its being an AGN, though it could be a variable 
background quasar. Since the SNABC compares all candi- 
dates to SN la and SN II-P spectral templates, it is, in effect, 
forcibly comparing all subtypes of CC SNe to SNe II-P. This 
leads us to believe that this SN is, in fact, a non-II-P CC SN. 
A similar situation is encountered for SNSDF0806.14. 



5.1.2 SNSDF0702.01 

SNSDF0702.01 was classified as a SN la [P(Ia) = 1], but 
with Xr = 13. At a separation of 3.61 ± 0.02 arcsec, this 
z = 0.18 transient is well offset from the centre of its spi- 
ral host galaxy, and so precludes the possibility of an AGN 
(though it could be a variable background quasar) . The high 
value arises from this object's R—i! colour, which does not 
fit the SN la template. As its absolute 7?-band magnitude 
is Ma = —17.01, we checked whether this could be a SN 
1991bg-like SN la by comparing its photometry to the N02 
SN 1991bg template. While the ^/-band magnitude matches 
the template, the R—i! and i' — s! colours do not. Though 
the z'-band magnitude and i! — z colour raise the possibil- 
ity that this is an early SN ITP, it is still too blue in the 
R band. We also checked whether the excess flux in the R 
band might be the result of a SN caught during shock break- 
out, by comparing the 7?-band photometry in our half-night 
stacks, but there was no discernible difference between the 
7?-band flux in the flrst two nights and in the second two. 
At this point we conclude that this object is too faint and 
too blue to be a SN la, and it might be either a very blue 
SN II-P, or a peculiar SN of a different kind. As detailed in 
Section [6] since this object is at z; = 0.18, it enters neither 
the SN la nor the CC SN rate calculations. 



5.1.3 SNSDF0702.30 

SNSDF0702.30 has two possible host galaxies, as shown in 
Fig. |6] a resolved galaxy designated hSDF0702.30a, and a 
compact galaxy to the NW (upp er right; hSDF070 2.30b). 



( Peng et al.||2010| ) to fit 



We used the software GALFII 

and subtract the larger galaxy, thus enabling us to perform 
photometry of each galaxy on its own. The resulting pho- 
tometry and best-fitting ZEBRA SEDs are shown in Fig. [T] 
Both galaxies agree well with the power-law SED of a star- 
forming galaxy at a high redshift (hSDF0702.30a at z = 2.0 
with = 3.5, and hSDF0702.30b at z = 1.7 with x^ = 0.8). 



http: / / users.obs.carnegiescience.edu/peng/work/galfit / galfit.html 
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Figure 7. Photometric redshift derivation for the two pos- 
sible hosts of SNSDF0702.30, shown in Fig. |6] The galaxy 
hSDF0702.30a appears on top, while hSDF0702.30b is below. 
Symbols as in Fig. |4] The first ^-PDF is peaked at z = 1.95, 
and the second z-PDF is peaked at z = 1.72. 



While the fit in Fig. [7] does not utilize UV data, the results 
agree with the nondetections observed in the FUV band, as 
seen in Fig. [6] 

Using the resulting z-PDF of hSDF0702.30a as a prior, 
the SNABC classifies this SN as a CC SN [P(Ia) = 0.04] 
at redshift z = 1.95, with xl = 37. The z-PDF of 
hSDF0702.30b, on the other hand, yields a different clas- 
sification: [P(Ia) = 0.68] at redshift z = 0.8, with Xr = 0.4. 
In this case, the SNABC chooses the smaller z-PDF peak 
at 2 ~ 0.8, instead of the main peak at z « 1.7, in or- 
der to avoid a high Xr value such as that achieved with the 
sharply peaked z-PDF of hSDF0702.30a. When run through 
the SNABC with a fiat z-PDF, the SN best resembles a 
CC SN [P(Ia) = 0.30] at z = 0.6 with x? = 1-0. The z-PDF 
constructed from the best-fitting redshifts of the other SN 
host galaxies does not change this result much; the posterior 
redshift changes to z = 0.7, with a lower Xr = 0.4. 

In this case, the SNABC is dominated by the SN II-P 
LF. Since the colours of the SN match those of a SN II-P, 
it places it at z < 1, the redshift range where the apparent 
magnitude of the SN would still match the SN II-P LF. 
This is also the reason it produces a high Xr value when 
forced to higher redshifts. In summary, SNSDF0702.30 may 
be either a CC SN a,t z = 0.6-0.8, or a non-la luminous SN at 
z = 1.7-1.95. The possible observation of overluminous non- 
la SNe at high redshifts in our sample is further discussed 
in Section |7.1.3| below. If it is a low-z CC SN, it will not 
be counted in the rates, as it is fainter than the detection 
limit adopted in Section [6] Since it may be a high-z non-la 
SN, we do not include this SN in our 1.5 < z < 2.0 SN la 
sample. 



5.1.4 SNSDF 07 05.17 

SNSDF0705.17 was classifled as a CC SN [P(Ia) = 0.02] at 
z = 2.87, with X? = 58. The offset of the candidate from 
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its host galaxy is 0.15 ± 0.10 arcsec, or ~ 1 ± 1 pixel. If 
one were to redshift a fiducial SN la template (i.e., at peak, 
with no stretch) to 2: = 2.87, the synthetic observed z'-hand 
magnitude would be z — 29.5, which is 3.9 mags fainter 
than the observed z — 25.6 ± 0.2 of the object. Thus, at 
this redshift, the object is too bright to be either a SN la or 
a normal CC SN. This object appears in epoch 4, which is 
separated from epoch 3 by only ~ 90 days in the observer's 
frame. In the object's rest-frame, this interval corresponds 
to ~ 23 days. The high redshift, coupled with the high Xr 
value, raises the suspicion that this candidate, even though 
it shows no variability in other epochs, is still an AGN. Al- 
ternatively, the object might be a hyperluminous SN Iln, or 
even a pair-instability SN. Since both luminous SNe Iln and 
pair-instability SNe decay slowly (e.g., Di Carlo et al.|2002 



with those of the z = 1.189 pulsational pair-instability SN 



Gal- Yam et al. 20091, if this object were one of the two it 



would likely have been detected in both epochs, unless it 
exploded between the two epochs. Our preferred conclusion 
is that this is an AGN and as such, we have removed it from 
our sample. 



5.1.5 SNSDF0705.18 

SNSDF0705.18 lies 3.06 ± 0.10 arcsec, about 3 half-light 
radii, from the closest (and only probable host) galaxy. 
We obtained a spectrum of this galaxy, which places it at 
z = 1.41. If this is indeed the SN's host galaxy, it is classi- 
fied as a SN la [P(Ia) = 0.98], with xl = 17. The SNABC 
is sensitive to the 2-PDF it receives as input, and since for 
this galaxy the input was a very narrow {a = 0.01) Gaussian 
centred on the measured spec-2:, we ran this SN through the 
SNABC once more, this time treating it as a hostless SN. 
This resulted in a classification as a CC SN [P(Ia) = 0.39] 
at a posterior redshift of 0.7, with a much better Xr = 0.2. 
At this redshift, the synthetic photometry derived from red- 
shifting the SN II-P template, at peak, would be 2' = 25.37 
mag. This is consistent with the measured z' = 25.6 ± 0.2 
mag. The z'-band master image of epoch 4 has a limiting 
magnitude of 27.24 mag. At z = 0.7, a galaxy would have 
to be fainter than —15.9 mag so as not to be detected. This 
could mean that the candidate is indeed a CC SN that went 
off in a dwarf galaxy undetected in the SDF (see, e.g., |Ar-| 
cavi et al.|2010[ ). Since the fit to a CC SN at z = 0.7 is much 
better than the earlier SN la classification, we treat this SN 
as a 'hostless,' intermediate redshift CC SN. As this SN is 
fainter than the detection limit adopted for this redshift bin 
(see Section[6] below), it will not be counted in the rates. To 
account for the possibility that this is a SN la in the range 
1.0 < 2 < 1.5, we add a systematic uncertainty of -f 1 to the 
number of SNe la in this bin. 



5.1.6 SNSDF0705. 30 and SNSDF0806. 35 

SNSDF0705.30 and SNSDF0806.35 are both classified as 
SNe la [P(Ia) = 0.90 and P(Ia) = 0.99, respectively] at 
high redshifts (2 = 1.93 and z = 1.94, respectively), but 
with high Xr values (34 and 22, respectively). While these 
SNe are both offset from the cores of their host galaxies 
(by 0.5 ± 0.1 and 0.7 ± 0.1 arcsec, respectively), they are 
much bluer than any of the SN la or CC SN spectral tem- 
plates. SNSDF0806.35 has R~i and i' — 2' colours consistent 



SCP 06F6 ( [Barbary et al.||2009| [Quimby et al.||2009[ ), red- 
shifted to 2 = 1.94. SNSDF0705.30, on the other hand, is 
even bluer. It might be a very blue non-la SN, or a back- 
ground variable quasar. As both of these SNe are clearly not 
SNe la, we exclude them from our 1.5 < 2 < 2.0 bin. 



6 DEBIASING: DERIVATION OF INTRINSIC 
SUPERNOVA TYPE AND REDSHIFT 
DISTRIBUTIONS 

The success rate of the SNABC depends on the intrinsic 
parameters of the SNe (e.g., type, age, redshift, and extinc- 
tion). P07a have found that degeneracies between these pa- 
rameters lead to misclassifications, which in this work may 
introduce biases in the SN rate calculations (i.e., if an ap- 
preciable number of SNe la are misclassified as CC SNe, the 
SN la rates will be systematically lower). In order to cor- 
rect for potential misclassifications, we follow the debiasing 
procedure described by P07a and P07b. We use the spectral 
templates from N02 to simulate a sample of 40,000 SN light 
curves, divided into four subtypes: la, II-P, Ib/c, and Iln. 
These templates have been normalised so that the _B-band 
absolute magnitude at maximum luminosity, for a stretch 



s = 1 (Perlmutter et al. 19991 SN la, is zero, in the Vega 



magnitude system. In order to construct the light curves in 
our sample, we follow the recipe outlined by [Sullivan et al.| 
( 2006 1 . For SNe la, the light curves are constructed accord- 



mg to: 



and 



m = 771^=0, s=i + Mb + M — ol[s — 1) 



ts 



(4) 



(5) 



where mz=o,a=i is the basic light curve, at 2 = and with 
s = 1, constructed from the spectral templates; Mb is the 
peak brightness in the B band, drawn from a Gaussian cen- 
tred on —19.37 mag, with a dispersion of ct = 0.17 mag, 
mimicking the intrinsic SN la dispersion in peak brightness 
( Hamuy et al.|1995| |1996| [Phillips et al.|1999[); /x i s the dis- 
tance modulus; a = 1.52 ±0.14 ( Astier et al.|2006 l; s is the 
stretch parameter of the SN, which is modeled as a Gaussian 
centred on s = 1 with a dispersion of a = 0.25, and allowed 



to vary in the range 0.7 5S s 5C 1.3 ( Sullivan et al.(2006) ; and 
ts is the age of the stretched-light-curve SN. 



The dispersion in s, taken from Sullivan et al. (20061, 



is larger than the observed dispersion among normal SNe la 
(e.g., [Howell et aL]|2007 l, in order to include both the very 
subluminous and overluminous SNe la. The above recipe re- 
sults in a LF that is consistent with those assumed by N06 
and 'Sullivan et al. (20061, and measured by Dilday et al. 
(20081). Recently, Li et al. (2011a I measured a larger frac- 



tion of subluminous SNe la than is represented here, which 
means our subsequent SN la rates may be underestimated. 



However, since the Li et al. (2011a I LF is not corrected for 



extinction, nor is it in a standard magnitude system, we can- 
not use it to estimate how many subluminous SNe la may 
be unaccounted for in our calculations. 

The CC SN light curves are constructed in much the 
same way, but without any stretching. Host extinction is 



added using the Cardelli, Clayton, & Mathis ( 1989 1 extinc- 
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tion law, with Rv = 3.1, and Av values drawn from the 
extinction model of N06: the positive side of a Gaussian 
centred on Ay = mag, with a dispersion of a = 0.62 for 
SNe la and a = 0.93 mag for CC SNe | | Sullivan et al.|2006 |. 

As with our observed SN sample, one sixth of the sim- 
ulated sample is assigned a random spec-z in the form of a 
Gaussian z-PDF with a = 0.01. The rest of the SNe in the 
sample are randomly assigned a redshift from the 2-PDF of 
the entire SDF, out to z = 3. Each simulated SN is assigned 
a 'real' redshift and a 'measured' redshift drawn from its 
z-PDF. This mimics the ZEBRA redshift determinations. 
While the simulated light curves are redshifted according to 
the real redshift, we keep the entire z-PDFs for the classifi- 
cation stage. 

The resulting light curves are 'observed' at a random 
day, and each measurement is assigned an uncertainty ac- 
cording to the photometric uncertainties measured in our 
survey. At redshifts z ^ 1, the light curves do not cover 
the full time period during which SNe could have been de- 
tected by the depth of our survey. One way to overcome 
this problem would be to extrapolate the light curves, but 
this might introduce systematic errors that are difficult to 
quantify. Instead, we chose to impose a flux limit on the SNe 
found in these bins; by raising the detection limit we nar- 
row the time period during which the SNe could have been 
observed, thus ensuring that we stay within the bounds of 
the observed light curves. 

In the 0.5 < z < 1.0 bin, the detection limit was raised 
to 25.0 mag in the z' band for all epochs. This reduces the 
number of SNe in this bin from 85 to 29, of which 26 are 
classified as SNe la and 3 as CC SNe. In the z < 0.5 bin, 
the necessary fiux limit leaves no SNe to work with; we thus 
cannot compute the SN rate in this bin. We note, however, 
that rates at z < 1 are much better measured by wider and 
shallower surveys that obtain light curves and spectroscopic 
confirmation for each SN (e.g., SDSS-II, SNLS). Our survey 
is designed specifically for detecting SNe at z > 1, and for 
classifying them with single-epoch photometry. 

We measure the success fractions of the SNABC in each 
epoch of observations by selecting only those SNe that would 
have been detected by our survey (i.e., those SNe which are 
brighter in the z' band than 26.3, 26.6, 26.4, and 26.7 mag 
for epochs 2 through 5, respectively, in the z > 1 bins, and 
brighter than 25.0 mag in all epochs for the 0.5 < z < 1.0 
bin), leaving 3,000 SNe from each subtype. The surviving 
SNe are then classified by the SNABC, and their redshift is 
determined as in Section [5] Next, the SNe are distributed 
into three redshift bins (0.5 < z < 1.0, 1.0 < z < 1.5, 
and 1.5 < z < 2.0), and the success fraction in each bin is 
calculated by dividing the number of correctly classified SNe 
by the total number of SNe in that bin. 

The resulting success fractions are used to calculate the 
probability of classifying a SN of any subtype as a SN la, as 
a function of the intrinsic distribution of SN subtypes (e.g., 
10 per cent SN la, 40 per cent SN II-P, 20 per cent SN Ib/c, 
and 30 per cent SN Iln). Using steps of 2.5 per cent, there 
are 12,341 possible distributions. In each redshift bin, and 
for each possible distribution, the SN la success fraction is 
computed by summing the fraction of SNe la that were clas- 
sified correctly, together with the fractions of CC SNe that 
were misclassified as SNe la. Each possible distribution is 
weighted according to the number of combinations in which 
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Figure 8. Observed (empty squares), flux-limited (filled circles), 
and debiased (solid line) SDF SN la and CC SN numbers. Filled 
squares denote the number of SNe in the 2 < 1 bins before appli- 
cation of the flux limit. SN la error bars arc la Poisson and classi- 
fication uncertainties, added in quadrature. CC SN 0.5 < z < 1.0 
debiased error bar is la Poisson and classification uncertainties, 
added in quadrature, and z > I debiased numbers are 2a upper 
limits (arrows). 



the different CC SN subtypes may be distributed for a given 
fraction of SNe la (i.e., if the fraction of SNe la is 50 per 
cent, there are many different combinations of CC SN frac- 
tions, whereas if the SN la fraction is 100 per cent, there is 
only one possible combination). 

After weighting the different distributions, we 
marginalise over all of the different combinations for a 
specific SN la fraction, and are left with the probability of 
classifying any SN as a SN la, as a function of the intrinsic 
SN subtype distribution. Using binomial statistics, this 
probability is used to answer the following question: Given 
the number of SNe classified by the SNABC as SNe la in a 
given redshift bin, the total number of SNe in that bin, and 
the probability of classifying any SN as a SN la, at a given 
intrinsic distribution, what is the most probable fraction of 
SNe la in our sample? From the resulting PDF we select the 
most probable value as the true fraction of SNe la in each 
redshift bin, and define the Icr uncertainty as the region 
that includes 68.3 per cent of the probability density. To 
this classification uncertainty we add, in quadrature, the 
statistical uncertainty, defined as the la Poisson uncertainty 
of the debiased number of SNe la in the redshift bin (or the 
Poisson uncertainty of the number of debiased CC SNe for 
the CC SN uncertainty). 

The raw and debiased distributions of SNe la and 
CC SNe are presented in Fig. [8] The debiased number of 
SNe la is the same as the raw number in the two z > 1 
bins, where our survey is mostly insensitive to CC SNe. The 
possibility that the z > 1 bins have been contaminated by 
luminous CC SNe (e.g., SNe Iln) is taken into account in 
the lower systematic uncertainty of the debiased number of 



SNe la in these bins: 28.0 



10.0 



-1-4.3,-1-0.0 
-3.1,-1.7 



+6.4,4 
5.3, 



at 1.0 < z < 1.5 and 



in the 1.5 < 2: < 2.0 bin. In the 0.5 < 2: < 1.0 



bin, the number of SNe la falls to 20. 3_ 



The post- 
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debiasing number of CC SNe, on the other hand, rises to 
8.7_3'5 _5'8. The errors for the above SN numbers are la 
Poisson and classification uncertainties, respectively. 



Lookback Time [Gyr] 



7 SUPERNOVA RATES 

In this section, we use the debiased distributions of SNe la 
and CC SNe to derive the SN la and CC SN rates in three 
redshifts bins: 0.5 < z < 1.0, 1.0 < z < 1.5, and 1.5 < 2 < 
2.0. Our rates are summarized in Table [3] and comparisons 
to the literature are given in Tables |4] and [5) and in Figs. [9] 
and llOl All rates from the literature have been converted to 
h = 0.7. In cases where they are originally reported in SNuB 
(SNe per century per IO^^Lq^b), we have converted them 
to volumetric rates using the redshift-dependent luminosity 
density function from BOS: 



jsiz) = (1.03 + 1.76 X 10* L0,B Mpc" 



7.1 The type la supernova rate 

The volumetric SN la rate is 



(6) 



flla({z)») 



(7) 



where {z)i is the effective redshift of each redshift bin i, Ma.i 
is the number of debiased SNe la in bin i, and tv{z) is the 
survey visibility time, integrated over the comoving survey 
volume element dV, at all redshifts z within bin i. 

The visibility time is the total amount of time we could 
have observed a SN, given the parameters of our survey. At 
a given redshift, we need to consider the dispersion in light 
curves that originates in three separate effects: the intrinsic 
dispersion in peak magnitude, the stretch-luminosity rela- 
tion, and the host-galaxy extinction. To account for these 
different effects, we calculate the visibility time of each pos- 
sible light curve, weight it by its probability (which is just 
the product of the probabilities of the separate effects), and 
sum over all possible combinations. 

As in the previous section, we construct each possi- 
ble light curve according to Equation |4] We construct light 
curves with all the possible combinations of peak magni- 
tude, stretch, and extinction, where Mb is allowed to vary 
as a Gaussian in the 2a range around —19.37 mag (where 
la = 0.17); the stretch parameter s is allowed to vary as a 
Gaussian centred on s = 1 with a dispersion of 0.25 in the 
range 0.7 ^ s ^ 1.3, with a = 1.52; and Av ranges between 
and 3 mag according to the N06 model. 

Each point in the light curve is multiplied by the ap- 
propriate detection efhciency taken from the functions in 
Section [3. 2[ and the entire light curve is then summed over 
the time it lies above the detection efficiency limit (the 50 
per cent detection efficiency limits for the z > 1 bins, and 
25.0 mag in the 0.5 < z < 1.0 bin). Finally, we sum over 
the different epochs (since for each epoch the detection effi- 
ciency limit is different), and end up with the visibility time 



1 2 3 4 5 6 7 

Scaled SFR 

o Various (see caption) 
V Neilietai (2007) 
Poznansl'i et ai. {2007b) 
- A Dalnien et ai. (2008) 



9 10 
Kuznetsova et ai. (2008) 
Rodney & Tonry (2010) 
SDF(R^=1.0) 
SDF (R„=3.1) 




0.8 1 1.2 

Redshift 



Figure 9. SN la rates from the SDF (filled squares) compared to 
rates from the literature. Circles denote \ow-z data from |Cappel- 
laro et al.| ||1999|l, [Hardin et al.| l|2000VPa in et al.| ||20 02l), 'Madg- 



et al. 


('2006'),' 


et al. 


(20081, 



wick et al.| l|2003|, [Tonry et al.| (2003), B lanc et al.[^2 004), Neill 
'C2006'),'Botticclla ct al.' (2 008( 1, |Dilday et al.[|[21308p , |Horesh 



2 008( 1, 
LOaTa 
alfpC 



Survey rates from jRo dney fc Tonry[ |[2010ll a re left-facing trian- 
gles. The GOODS rates from [Dahlen et aT] l [2008[ l are denoted 
by upturned triangles. Right-facing t riangles are for |Kuznetsova| 
[et al.[ | [2008l l. Our initial SDF results ( [Poznanski et al.|2007b^ are 
shown as diamonds. Filled squares (circles) denote the SDF rates, 
derived with an extinction law with Ry = 3.1 {Ry = l).The cos- 
mic SFH from Y08, has been scaled to fit the low-2 data. The 
shaded area denotes the plausible range of SFHs with power-law 
slopes between 3 and 4, out to ^ = 1, and between —2 and for 
z > I. All vertical error bars include statistical and systematic 
uncertainties added in quadrature. Horizontal error bars indicate 
the SDF redshift bins. 



of our entire survey. Symbolically, 

U{z) ^ ^ JJJ dhlB ds dAv p{Mb) p{s) p(A„) 

epoch 



X / €[mz(t)]dt. 



(8) 



We take the weighted average of the redshifts in a bin as 
the bin's effective redshift, where the weight is the visibility 
time integrated over the volume element within that bin: 



jU z dV 
jU dV ' 



(9) 



The uncertainties of the rates are the classification and 
Poisson uncertainties of the debiased SN la numbers, prop- 
agated and added in quadrature. To test how the uncer- 
tainties in the detection efficiency functions, as plotted in 
Fig. [2j affect the rates, we ran 500 Monte Carlo simulations 
in which each efficiency measurement was varied according 
to its uncertainty. This produced variations in the detec- 
tion efficiency limits of ±0.1 mag. This propagates to a la 
dispersion in the SN la rates that is between one and two or- 
ders of magnitude smaller than our main uncertainties, thus 
having a negligible effect on the resulting rates. The SN la 
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Table 3. SN la and CC SN numbers and rates 



Subsample 


0.0 < z < 0.5 


0.5 < z < 1.0 


1.0 < z < 1.5 


1.5 < 2 < 2.0 


Total 


25 


85 


28 


12 


SN host galaxies with spec-z 
Hostless SNe 


4 



13 
11" 


7 

1 






SNe la (raw) 

SNe la (after flux limit) 

SNe la {debiased)'^ 

SN la rate (10""' yr"! Mpc'^) 

SN la rate without host-galaxy extinction 

Effective redshift 


7 



64 
26 

ori Q-t-5.2,+5.8 
^U.O_4 7 g 3 

"•'^-0.41 

60"'"'''2^ 
"■""-0.31 

0.74 


28 
28 

^^'■U-S. 3,-7.6 
n f,A+0.2h 
U.84_p 28 

62+°" 
"■"■^-0.21 

1.23 


lO*" 
10 

iQ (,-^4.3, +0.0 

^•^^-0.37 
45+0-20 

1.69 


CC SNe (raw) 

CC SNe (after flux limit) 

CC SNe (debiased)<* 

CC SN rate (lO--* yr"! Mpc'^) 

CC SN rate without host-galaxy extinction 

Effective redshift 


18 



21 
3 

o „+3.2,+9.3 
»-'_3.5,-5.8 

0.66 





< 3.8,-1-20.2 





< 3.8,-1-4.7 



"This includes SNSDF0705.18, which is treated as a hostless SN, as detailed in Section 5.1.5 

''Two of the 12 SNe in this bin are clear non-la transients. 

'^Errors are Icr Poisson and classification uncertainties, respectively. 

'^Errors in the 0.5 < z < 1.0 bin are la Poisson and classification uncertainties, respectively. 

The z > 1 rates are upper limits. Errors are 2a Poisson and classification uncertainties, respectively. 



rates, with and without taking host-galaxy extinction into 
account, are shown in Table |3] 



7.1.1 High-redshift dust 

As star formation increases with redshift, so does injection 
of dust into the interstellar medium, leading to an expected 
increase of extinction with redshift (e.g., [Mannucci, Delia] 
Valle, fc Panagia|2007 1 . This effect should lead to a decrease 
in the number of observed SNe at high redshifts. [Mannuca] 
|et al.| ( [2C)07[ ) have shown that at high redshifts {z > 1) a large 
fraction of SNe, both CC SNe and SNe la, would be missed 
in optical surveys, due to extinction by dust in massive star- 
burst galaxies, which make up a larger fraction of the galaxy 
population at higher redshifts \Le Floc'h et al.|2005 Daddi 



et al. 



et al. 



2005 Perez-Gonzalez et al.|2005 I . Using the Mannucci 
( |2006| ) DTD model, [Mannucci et aL] ^Q7\ calculated 
that in the range 1.0 < z < 2.0, 15 to 35 per cent of SNe la 
would be missed. Assuming a power-law DTD model with 
a slope of —1 (see Section [s] below), the fraction of missing 
SNe la would be 5-13 per cent in the above redshift range 
(F. Mannucci, private communication). The corrected SN la 
rates are shown in Fig. [9] and in Table [4] 



7.1.2 Different extinction laws 

Throughout our classification, debiasing, and subsequent 
derivation of the SN la rates, we have assumed a |Cardelli] 
et al. ( 1989 1 extinction law with the Galactic average Rv = 



3.1. However, several SN surveys have discovered SNe la 
that underwent extinction best fit with lower values of Rv, 
from 1.7 to 2.5 ( |Tripp |1998}|Hicken et al.|2009[[Wang et al 



2009). To gauge the systematic effect of lower Rv values on 



our rates, we reran the classification, debiasing, and rates 
derivation stages assuming an extinction law with Rv = 1. 



The resultant rates are shown as filled diamonds in Fig. [9] 
They are consistent with the rates derived with Rv =3.1, 
but in the three redshift bins are systematically lower by 6, 
26, and 43 per cent, respectively. 

Whereas it is possible that the extinction law in the 
immediate vicinity of SNe la is different from the Galactic 
average, recent studies (e.g., |Guy et al.|2010||Foley fc Kasen| 
20111 have raised the possibility that SNe la have an in- 



trinsic colour scatter, which together with dust extinction is 
responsible for their overall reddening. Chotard et al. (2011 1 



have found that once they corrected for an intrinsic scatter 
in the Si and Ca features of 76 SNe la spectra, they recovered 
a Cardelli et al. ( 1989 1 extinction law with Rv ~ 2.8 ± 0.3, 



consistent with the Galactic average value. We do not add 
the systematic uncertainty produced by different values of 
Rv to our final quoted SN la rates. However, in Section [8] 
we do take this systematic uncertainty into account when 
deriving the SN la DTD. 



7.1.3 Contamination from high-z non-la transients 

While our survey is largely insensitive to CC SNe ai z > 
1, there remains the possibility of contamination by non-la 
luminous SNe (e.g., [Smith et al.||2007| [Quimby et al.|[2007 



Barbary et al. 2009[ ). As described in Section 5.1 we have 
discovered two luminous non-la SNe in the 1.5 < z < 2.0 
bin. This ratio of 2:10 non-la SNe to SNe la is consistent 
with the 1:11 ratio found by Barbary et al. (20101, who 



found one non-la transient (SCP 06F6) and 11 field SNe la 
in the redshift range z = 0.6-1.3. 

As for contamination by AGNs, the extremely 
blue colours of SNSDF0705.30 and the classification of 
SNSDF0705.17 as a CC SN at z = 2.87 hint that these 
objects might be variable background quasars, as discussed 
in Section [5.1 1 This is consistent with the expected number 
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Table 4. SN la rate measurements 



Redshift 



A^la Rate [10-* yr-i Mpc" 



Reference 



0.01 
< 0.019" 
0.0375 

0.09 
0.098 

0.1 

0.13 

0.14 

0.15 

0.15 

0.2 

0.2 

0.25 

0.25 

0.3 

0.3 

0.35 

0.35 
0.368 

0.40 

0.45 

0.45 

0.46 
0.467 

0.47 



70 

274 
516= 

17 

19 
516= 

14 

4 
516= 
1.95 

17 
516= 
1 

516= 
31.05'' 
516= 

5 
4.01 

17 
5.44 

9 
5.11 

8 

73 



0.183 ±0.046 

n acc+0. 034, +0.043 
U.ZDO_Q Q 

n 978+0. 112, +0.015 
"■^'"-0.083,-0.000 



0.293 
0.24j 



-,+0.09 
="-0.07 
1+0.12 
-0.12 
n 9c;q+0. 052, +0.018 
U-^Ja-0. 044, -0.001 
n 1 co+0. 056, +0.035 
U.10S_Q Q43 Q 035 

„ „„+0.22,+0.07 
U.^O_0.i3 -0.04 
n on7+0-''38, +0.035 

U.OU/_g Q34 Q 

p. or,+0.23,+0.07 
"••-'^-0.23,-0.06 
n 1OQ+0.042, +0.018 I n 49 
U.189_Q P34 ± 0.42 

n 0/10+0.032, +0.082 

"■•340_Q Q3Q _Q 

0.17 ± 0.17 

n occ+0. 031, +0.182 
"■•jDO_g Q g]^2 

n 0^+0.16, +0.21 
U-J^-0.15,-0.22 
„ .„.+0.03'?, +0.396 
"■^■J^-0. 034, -0.016 

0.530 ± 0.024 

n 0^+0.19, +0.07 
"••""^-0. 19, -0.03 
„ „-. +0.05, +0.08 
"••^J^-0.05,-0.03 

n 50+0.39 

"•"•^-0.17 

0.73 ±0.24 

I +0.15, +0.12 
-0.15,-0.04 

0.48 ±0.17 

n -9+0. 06, +0.13 
"•^^-0.06,-0.09 
-1+0. 37, + 1.66 
-0.27,-0.26 



0.31^ 



0.80^ 



Cappellaro et al. 1 1999 1 



Li et al. 



(2011b I 



Dilday et al. 



Dilday et al. 



I 2010a I 



I 2008 I 



Madgwick et al. 



Dilday et al. 



Blanc et al. 



Hardin et al 



Dilday et al. 



(20031 



I 2010a I 



(20041 



(20001 



I 2010a I 



Rodney &; Tonry 



Horesh et al. 



Dilday et al. 



(20101 



( 2008 I 



I 2010a I 



Dilday et al. 



Barris fc 'iSnry' 1^006 i 



I 2010a I 



Botticella et al. 



I 2008 1 



Dilday et aT] |26ioa| 



Bam s Ic Tonry ( 2006 1 
Rodney fc Tonry^ ( ^201o| 



Neill et al. l|2007ll 

Kuz netsova ct al. {'2008') 
Barris & Tonry ( 2006 1 
Rodney & Tonry ( 20To| 
~5nry et al . ^ (,2003| ) 



Neill et al. 



I 2006 1 



Dahlen et al. 



I 2008 I 



0.55 
0.55 
0.55 
0.552 
0.65 
0.65 
0.714 
0.74 
0.74 
0.75 
0.75 
0.80 
0.83 
0.85 
0.95 



38 
29 

6.49 

41 

23 
10.09 

42 

5.5 
20.3 

28 
14.29 
18.33 

25 
15.43 
13.21 



0.568 



+0.098, +0.098 



0.088,-0.088 

2.04 ± 0.38 

n 09+0. 14, +0.07 
"••-'^-0.14,-0.07 
+0.10;+0.26 
0.10,-0.27 

1.49 ±0.31 

+0.17, +0.14 
0.17,-0.08 
+0.19, +0.54 
0.19,-0.70 
+0.36 
0.32 
+ 0.33 
0.41 

1.78 ±0.34 

«o+0.21,+0,23 
^°-0.21,-0.14 



0.63_ 
1.4 

0.49^ 

1.13j 
0.433 
0.79"* 



Pain ct al.'('2002')" 


Barris & Tonry ( 


2006jl 


Rodney & Tonry 


(20101 


Neill et al.'l'2007l) 


Barris & Tonry ( 
Rodney & Tonry 


2006) 

(20101 


Neill et al. i 2007 




Poznanski et al. 


2007b 







n 0-5+0.25 
u.»o_g 25 



1.30 
0.78 
0.76 



+0.33, +0.73 
0.27,-0.51 
+0.22, +0.31 
0.22,-0.16 
+0.25, +0.32 
-0.25,-0.26 



Barris & Tonry ( 2006 1 


Rodney & Tonry 


(2010') 


Kuznetsova et al. 


{,2008J 


Dahlen et al. 


(20081 


Rodney &; Tonry 


(20101 


Rodney & Tonry 


(20101 


Rodney & Tonry 


(2010) 


Kuznetsova et al. 


12008 


Dahlen et al. 


2008|l 


Poznanski et al. I 2007b 



1.05 
1.20 
1.21 
1.23 
1.23 



11.01 
8.87 
20 
10.0 
28.0 



0.28, +0.36 
0.28,-0.41 
+0.35 
0.30 
+0.36, +0.38 
0.29,-0.32 
+0.45 
0.56 
0.84^-- 



0.79 

0.75 
1.32 

1.05 



SDF (this workj 



1.55 
1.61 
1.67 
1.69 



0.35 
3 

3.0 
10.0 



0.12 
0.42 
0.81 
1.02 



+0.58 
0.12 
+0.39, +0.19 
0.23,-0.14 
+0.79 
0.60 
+ 0.54 
0.37 



Kuznetsova et al.' (12008' 


Dahlen ct al. 


(2008) 


Poznanski et 


al. 




2007b 



SDF (this workJ 



Note — Redshifts are means over the redshift intervals probed by each survey. A'^i^ is the number 
of SNe la used to derive the rate. Where necessary, rates have been converted to h = 0.7. 
Where reported, the statistical errors are followed by systematic errors, and separated by commas. 
The uncertainties of the SDF results are statistical and systematic, added in quadrature. 
(|2011b| consider SNe la within 80 Mpc. 



" Li et al. 

Rates have been con verted to volumetric rates using Equation 1 
Dilday et al. (2010a) compute their rates using 516 SNe la in tne redshift range z < 0.5. 



Botticella et aTi (|2008|) found a total of 86 SN candidates of all types. See their section 5.2 
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Table 5. CC SN rate measurements 



Redshift 



Ncc Rate [10^* yr^i Mpc"^] Reference 



< 0.0066" 
0.01 

< 0.014" 
0.21 
0.26 
0.3 
0.3 



92 
67 
440 
44.95^= 
31.2'' 
17 
117 



> 0.96 
0.43 ±0.17 
0.62 
1.15 

1.88 
2.51 
1.63 



+0.07, +0.17 
0.07,-0.15 
+0.43, +0.42 
0.33,-0.36 
+0.71 
0.58 
+0.88, +0.75 
0.75,-1.86 
+0.34, +0.37 
0.34,-0.28 



Smartt et a l. (2009) 
Cappellaro et al.^ ( ^1999| '' 



Li et al. 



(2011b 



Dahlen et al 



Bazin et al. 



1, et al.' ("200811 
:o et al ^0(j5| l'' 

falTproli 



Botticella et al. 
Cappellaro ct aL^_^l 



( 2009 I 



0.66 

0.7 



8.7 

17 



6.9 



-9.9 
-5.4 



3.96 



+1.03, + 1.92 
-1.06,-2.60 



SDF (this work) 
Dahlen et al.|l|2004l 



Note - Where reported, the statistical errors are followed by systematic errors, and separated by commas. 
The uncertaint ies of the S DF resul t s are st atistical and systematic, added in quadrature. 
(2009) and 



' Smartt et al. 



same as m 



febleli] 



Li et al. 



(2011b I consider CC SNe within 28 and 60 Mpc, respectively. 



^Total number of CCT SNe and SNe la detected throughout the survey. 



of contaminating AGNs in our sample, as detailed in Sec- 
tion |4.1[ In summary, beyond the non-la objects we have 
identified, contamination of the 1.5 < z < 2.0 SN la sample 
is unlikely. 



7.1.4 Biase,s in the photometric redshifts 

As shown in Fig. [3] there is a bias in our photo-z method 
towards overestimation of the redshift at high redshifts. 
This is caused by the dearth of spectroscopic redshifts at 
1.5 < z < 2.0 (only 24, or ~ 6 per cent, of the training-set 
galaxies), and the inherent difficulty of measuring the red- 
shift of late-type galaxies (see Section 4.2). We have taken 



spectrum, introducing a possible systematic uncertainty into 
any survey (such as the current work) which probes this part 



two steps to overcome this bias. First, the measured colours 
of the SNe were compared to those predicted by the tem- 
plates of SNe la, SNe II-P, SNe Ib/c, and SNe Iln, at differ- 
ent redshifts, spanning the entire < 2 < 2 range. 

Eight out of the ten 1.5 < z < 2.0 SNe 
agree within 2a with the fiducial SN la template col- 
ors one would observe at their host galaxies' photo-z 
(namely, SNSDF0705.21, SNSDF0806.31, SNSDF0806.46, 
SNSDF0806.50, .25, SNSDF0705.29, SNSDF0806.38, and 
SNSDF0806.57). SNSDF0702.28 may be an example of the 
bias seen in Fig. [3] whereas its late-type host galaxy has 
a photo-z of ~ 2, the SN colours favour a lower redshift 
of ~ 1.6-1.7. Finally, the colours of SNSDF0806.32 favour 
the SN Iln template over the entire 1.2 < 2: < 2.0 redshift 
range. The possibility that this SN has been misclassified as 
a SN la is taken into account in the systematic uncertainty 
of the SN la rate in this bin, as quoted in Table [3] 

To further investigate the issue of the redshifts of the 
candidate z > 1.5 SNe, and to check whether any of them 
are contaminating AGNs, we are pursuing HST and ground- 
based spectroscopic observations of these hosts. 



7.1.5 Probmg the UV part of the SN spectrum 

From a theoretical standpoint, the spectra of SNe la at high 
redshifts may differ from their low-redshift counterparts due 
to changes in, for example, progenitor metallicity. Such dif- 
ferences are expected to show up in the UV part of the SN la 



of the spectrum (Hoeflich, Wheeler, & Thielemann 1998 



Lentz et al.||2000 Sauer et al.||2008 1. Several recent surveys 
have found evidence for such differences between low- and 



high-redshift SNe la (e.g., Kessler et al.|2009 Cooper et al. 
2009 Foley et al.||2010" ), which might provide an additional 
explanation for the high values of the two peculiar SNe 
in our 1.5 < z < 2.0 sample. 



7.2 The core-collapse supernova rate 

Since our survey is insensitive to normal CC SNe at redshifts 
higher than 1, we do not use the debiased results to derive 
the rates in the 1.0 < z < 1.5 and 1.5 < z < 2.0 redshift 
bins. We now proceed to derive the CC SN rate in the 0.5 < 
2 < 1.0 redshift bin. 

To account for the division of CC SNe into subtypes, 
in the calculation of the visibility time we have weighted 
the contribution of each subtype according to its fraction of 
the total CC SN population, and then summed the differ- 
ent contributions. The CC SN subtype fractions were taken 
from the volume-limited sample of Li et al. (2011a I, with 



two alterations: (a) the SN II-P and SN II-L fractions have 
been combined, as the separation between these subclasses 
is currently ill-defined ( Poznanski et al^|2002 l; and (b) the 
SN Ib/c and SN lib fractions have also been combined, since 
their light curves are nearly identical ( Benson et al.|[T994 |. 
The final volume-limited CC SN fractions are 60.0 per cent 
II-P/L, 33.5 per cent Ib/c/IIb, and 6.5 per cent Iln. We note 
that Li et al. ( 2011a|b l only targeted ~ L* galaxies, and so 
the CC SN fractions and rates might be different for an un- 
targeted survey (e.g., Arcavi et al.||20i0 l. We calculate the 
flux-limited fractions at the effective redshift of 2 = 0.66 as 
being 37 per cent II-P/L, 44 per cent Ib/c/IIb, and 19 per 
cent Iln. 

As in the previous section, the visibility time of each 
CC SN subtype was derived using Equation |8] but without 
stretch. In the present case, Mb was limited to the 2cr range 
around the peak magnitude of each subtype. The probability 
for Av was drawn from a one-sided Gaussian PDF centred 
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Table 6. SN rate uncertainty percentages 



B 
re 
cc 

z 

CO 

o 
o 



Scaled SFR 






Cappellaro etal. (1999) 




a Dahlen et al. (2004) 




- □ Cappellaro et al. (2005) 




Bottlcella et al. (2008) 




V Bazin et al. (2009) 




o Smartt et al. (2009) 




< LI etal. (2011a) 






- ♦ SDF(R^=1.0) 






. SDF(R^=3.1) 















0.4 0.6 

Redshift 



Figure 10. CC SN rate from the SDF (filled square) compared 
to rates from the literature: [Cappellaro et al.||1999| (right-facing 
triangle) , [Dahlen et al.|2004| (upward triangles), [Cappellaro et al. 
[2005[ (square), [Botticella et al.|2008[ ( diamond), [Bazin et al.|2009 
(downward triangle), lower limit from [Smartt et al.|2009[ (circle), 
and [Li et al.[|2011b| (left-facing triangle). As in Fig. [9| the SFH 
from Y08 has been scaled to fit the low-2: data. All vertical error 
bars from the literature are la uncertainties. The horizontal error 
bar indicates the SDF redshift bin. 



on with a dispersion of a = 0.93, and the probability 
for Mb was drawn from the LF of each subtype. Without 
host-galaxy extinction, the rates of each CC SN subtype (in 
units of 10"^ SNe yr'^ Mpc"^) are 1.3 for SNe II-P/L, 0.4 
for SNe Ib/c/IIb, and 0.1 for SNe Iln. This resuhs in an 
overall rate of Rcci{z) = 0.66) = 1.8^11 x lO"* SNe yr'^ 
Mpc~'^. Once host-galaxy extinction is added, the rates of 
each CC SN subtype (in the same units) become 5.8 for 
SNe II-P/L, 0.9 for SNe Ib/c/IIb, and 0.2 for SNe Iln. This 
yields an overaU rate of Rcc{{z} = 0.66) = 6.9l5;4 x 10"" 
SNe yr"^ Mpc"^. After correcting for the fraction of CC SNe 
missed due to high-redshift dust ( [Mannucci et al.|2007[ ), the 
final CC SN rate is Rcc{{z} = 0.66) = 6.9+^;^ x lO"'^ SNe 
yr~^ Mpc"''. This value is consistent with both the rate 
derived by D04 in this redshift bin and with the scaled Y08 
SFH at that redshift, as shown in Fig. |10[ We present a 
summary of CC SN rates from the literature, along with 
our measured rate at {z) — 0.66, in Table [s] 

The statistical and systematic uncertainties affecting 
the SN la and CC SN rates are summarised in Table |6l 



8 THE TYPE lA SUPERNOVA DELAY- TIME 
DISTRIBUTION 

In this section we make use of our measured SN la rates, to- 
gether with published rates at various redshifts, to recover 
the DTD. The different SN la rates used in our fits are 
presented in Table [4| W here necessary (Cappellaro et al 



1999 



2003 



IHardin et al.[[2000[ [Pain et al.[[2002: Madgwick et al 
[Blanc et al.|2004 1 , rates from the literature have been 



converted to volumetric rates using the redshift-dependent 
luminosity density function from BOS (see Equation [6|. Fur- 
thermore, all rates have been converted to h — 0.7. We make 



Uncertainty 


0.5 < z 


< 1.0 


1.0 < z 


< 1.5 


1.5 < z 


< 2.0 


SN la rates 


Poisson 


+25/ 


- 23 


+23/ 


- 19 


+42/ 


- 31 


Classification 


+28/ 


- 45 


+3/- 


- 27 


+0/- 


- 18 


High-2 dust 


+3/ 


- 


+6/ 


- 


+9/ 


- 


Extinction law" 


+0/ 


- 6 


+0/- 


- 26 


+0/- 


- 43 


Total 


+41/ 


- 51 


+29/ 


- 33 


+51/ 


- 36 



CC SN rates 



Poisson 


+37/ - 40 


Classification 


+107/ - 67 


High-2 dust 


+32/ - 


Extinction law" 


+0/ - 52 


Total 


+145/ - 78 



All uncertainties are reported as percentages of the rates. 
"This uncertainty is not added to the final quoted rates, but is 
propagated direct ly into the derivation of the SN la DTD 
(see Section 7.1.21. 



use of all the SN la rate measurements in Table[4] except for 
the Barris & Tonry ( 2006 1 measurements, which have been 
superseded by [Rodney fc Tonry[ ( [20Tol ); the [Kuznetsova et al.[ 
(2008) measurements, which make use of much the same 
data as D08; and our initial results reported by P07b, which 
are superseded by the present results. In total, there are 36 
SN la rate measurements, of which 31 are at z < 1 and 5 at 
z>l. 

We recover the DTD by convolving different trial DTDs 
with various SFH fits from the literature, resulting in a 
model SN la rate evolution. One such SFH is the one pre- 
sented in fig. 2 of HB06 (HB06c). Other recent estimates 
of the SFH and their parametrizations (e.g., Y08 and O08) 
can be approximated by broken power laws, with a break 
at 2 — 1, and various power-law indices above and be- 
low the break. To test the systematic uncertainty in our 
DTD derivation produced by this range of possible SFHs, we 
parametrize the SFH as being proportional to (l + z)^, with 
7 in the range 3-4 at z < 1, a break at z — 1, and 7 values 
in the range —2 to at z > 1. This range of parametriza- 
tions covers most of the SFHs that have been recently pro- 
posed. The indices, breaks, and normalizations of each SFH 
at z = are collected in Table [7j and the resulting SFHs 
are shown in Fig. [11| For a given SFH, variations of the 
normalization will translate to inverse scalings of the am- 
plitude of the best-fitting DTD, without affecting the DTD 
shape, which is our main interest here. There remains con- 
siderable debate among different authors as to the amount 
and the redshift dependence of extinction corrections in SFH 
estimates (see, e.g., [Bouwens et al.[|2010[ [Robertson et alT] 
[201C)j ) . Different extinction correction choices can shift much 
or all of a SFH curve up or down by up to a factor of two 
(F. Mannucci, private communication). To account for this 
uncertainty, we also calculate the range in DTD amplitude 
that results when the SFH varies between the extreme case 
of O08u and the HB06c level. 

Throughout this derivation we assume a 'diet-Salpeter' 
initial mass function (IMF; [BeU et al.[[2003[ ). This IMF as- 
sumption means that the SFHs of HB06 and Y08, who as- 
sumed a Salpeter (19551 IMF, are scaled down by a factor 
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Figure 11. SFH measurements and parametrizations. Measure- 
ments include tiie compilation from HB06 (circles) and addi- 
tional data presented by [Horiuchi &: Beacom| | |2010 i (squares; 
here noted as HBIO). The dashed line represents the Cole et al.| 
l |2001[ l parametrization with parameter values from HB06. The 
solid (Y08), dot-dashed (O08u), and dotted (O081) lines are power 
laws with parameter values as detailed in Table [7] 
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Figure 12. Best-fitting values and 68 per cent statistical un- 
certainties of the slope /3 of a power-law DTD of the form 
<!/(<) = \E'i(t/l Gyr)'', when convolved with various SFHs, as 
marked. See Table [7] for SFH abbreviations and parameters. 



Figure 13. Top panel: Observed SN la rates compared to predic- 
tion from convolution of the Y08 SFH with a best-fitting power- 
law DTD of the form ^(t) = 1'i(Vl Gyr)'' (solid line). Non- 
independent measurements, which are therefore excluded from 
the fits, are not plotted — [Kuznetsova et al.| | |2008[ l and P07b, 
which are superseded by D08 and this work, respectively. The 
shaded area is the confidence region resulting from the 68 per 
cent statistical uncertainty of /3 from the convolution of the DTD 
with the Y08 SFH fit. Bottom panel; Same as top panel, but for 
each of the SFHs in Table [7] and showing the combined effect of 
the 68 per cent statistical uncertainties of /3. The thin dashed lines 
indicate the 68 per cent uncertainty region produced without the 
new SDF measurements. 



of 0.7. We then use the x statistic to find the best-fitting 
values of the parameters of the DTD, along with their sta- 
tistical 68 and 95 per cent confidence regions, defined as the 
projections of the Ax^ ~ 1 contour on the parameter axes. 
Below we present reduced values, denoted by Xr- To the 
statistical uncertainty of the parameters we add the system- 
atic uncertainty that originates in the shapes of the differ- 
ent SFHs. Finally, for each model we calculate the number 
of SNe la per formed stellar mass, integrated over a Hubble 
time. 

We first test a power-law DTD of the form — 
'i>i{t/l Gyr)''. Such a power law, with 13 ~ —1, is generic 
to the DD scenario, where two WDs merge due to loss of 



energy and angular momentum to gravitational waves (see, 
e.g., Maoz et al.|2010 l. Several recent experiments, in differ- 
ent environments and different redshifts, have indeed found 



best-fitting DTDs consistent with this fo rm (Totani et al. 
Maoz fc Badenes||2010] |Maoz et al.]|2010[ |2011[ ). The 



2008 



DTD is set to zero for the first 40 Myr, until the forma- 
tion of the first WDs. We fit for the normalization ^'i and 
the slope /3. Based on the Y08 SFH fit, we find best-fitting 
values of P — 1.1 ± 0.1(0.2), where the statistical uncer- 
tainties are the 68 and 95 (in parentheses) per cent confi- 
dence regions, respectively. The range of SFHs tested here 
adds a systematic uncertainty of Ig lo- The best-fitting val- 
ues of P for all four SFHs, with their respective reduced 
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values, appear in Table [t] These best-fitting values re- 
sult in reduced Xr values of 0.7 to 0.8, for 34 degrees of 
freedom (DOF) for all SFH fits. The number of SNe la per 
formed stellar mass, integrated over a Hubble time, lies in 
the range A^sn/M* = (0.5-1.5) x 10"^ M0~\ where this 
range accounts for the statistical uncertainties in P and 'I'l. 
However, the uncertainty in the normalization of the SFH 
is such that this range might easily be higher by a factor of 
two (F. Mannucci, private communication). The best-fitting 
values of j3 are presented in Fig. |12| and the resulting SN la 
rate evolution tracks are presented in Fig. |13[ 

Whereas the power law discussed above extends all the 
way back to t = 40 Myr, it is possible that at early times the 
DTD is dictated not by the WD merger rate, but rather by 
the supply of progenitor systems. Pritchet, Howell, fc Sul-| 



Lookback Time [Gyr] 

01234567 8 9 10 11 



livan 



(20081 have suggested a t power-law DTD, which 
is proportional to the formation rate of WDs. A pure t~^^^ 
power law, convolved with the HB06c, Y08, and O081 SFHs, 
produces fits with a minimal Xr > 1-5 for 35 DOF, ruling 
out this model at the 95 per cent confidence level. The O08u 
SFH results in a fit with a minimal Xr value of 1.4, which 



is marginally acceptable. Matteucci et al. (20091 also argue 



against this model, as it does not reproduce the observed 
G-dwarf metallicity distribution in the solar vicinity (see 
their fig. 7). The resulting SN la rate evolution tracks are 
presented in the top panel of Fig. |14[ 

Another possibility is that the DTD is controlled by 
the WD formation rate up to some characteristic time tc, 
beyond which newly formed WDs no longer have the com- 
bined mass to constitute SN la progenitors; from this point 



on only the merger rate sets the DTD. The Greggio (20051 
DD3-close model, for example, is such a broken power law, 
with t~^^^, t~^''^ , and a break at tc — 0.4 Gyr. This value for 
tc corresponds to the lifetime of 3 Mq stars. A larger value 
of tc would imply that WD binaries with a smaller primary 
mass can explode as SNe la, and ultimately contribute to the 
observed SN la rate. We therefore investigate whether the 
SN la rate data may be fit by a broken power law behaving 
as ai t < tc, and as t~^ thereafter. Fitting for tc and 

the normalization ^Pi, we find that tc lies in the 68 per cent 
confidence range 0.04-0.48 Gyr. As a t~^^'^ power-law DTD 
at all times is still an acceptable option for the O08u SFH, 
we cannot constrain tc at the 95 per cent confidence level 
for that SFH. However, the other SFHs suggest that tc may 
be lower than ~ 0.8 Gyr, at the 95 per cent confidence level. 
The best-fitting parameters result in reduced Xr values of 
0.7-0.9, for 34 DOF for all SFH fits. The integrated number 
of SNe la per stellar mass formed resulting from this DTD 
lies in the range Nsn/M* = (0.5-1.0) x 10"^ Mq"\ where 
the uncertainty derives from the normalizations of the SFHs 
and from the statistical uncertainty tc- This range is simi- 
lar to that obtained with the single power-law DTD. The 
best-fitting parameters, along with reduced x^ values, are 
presented in Table [7| and the resulting SN la rate evolution 
tracks are presented in the centre panel of Fig. | 14[ 

Finally D04, D08, and|Strolger et al.| p0"04[ |2010[) ad- 



vocate a Gaussian DTD with parameters r = 3.4 Gyr and 
a = 0.2r. D04 used the SFH determined by |Giavalisco et aL] 
(20041 in order to derive the parameters of the Gaussian 



DTD. As we use different SFHs, we leave the normaliza- 
tion of the Gaussian, ^g, as a free parameter. The best- 
fitting value, derived with the HB06c SFH fit, has a mini- 
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Figure 14. Observed SN la rates compared to predictions from 
convolution of the SFHs in Table[7]with a best-fitting (top) power- 
law DTD of the form *(t) = <I'i(t/l Gyr)-^^^; (centre) broken 
power-law DTD of the form *(t) oc J-^/^ up to tc, and >!>(<) oc 
afterward; and (bottom) DOS Gaussian DTD. Symbols are as 
marked. 
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Table 7. Star formation histories and resultant best-fitting DTD parameters. 





SFH 


Power-law DTD 


Broken power-law DTD" 


Ref.'' 


Parametrization'^ 


xVdof 


tc [Gyr]= x7DOF 


Galactic dust extinction: Ry = 3.1 



HB06c 

Y08 S(0) 

O081 S(0) 

O08u 5(0) 



Cole et al. 



1 2001 1 with values from HB06 
3.4, Zi = 1, 72 = -0.3 

3, ^i, = 1, 72 = -2 

4, zi, = l, 72 = 



1.11 



-0.10(0.24) 



17.8, 71 
17.8, 71 
17.8, 71 



-0.10(0.20) 

1.1 ± 0.1(0.2) 



1.23 
0.96 



+0.16(0.40) 
-0.13(0.25) 
+0.09(0.17) 
-0.07(0.16) 



0.7 
0.7 
0.8 
0.7 



+0.14(0.70) 
"■"'-'-0.01(0.01) 
„ „.+0.17(0.72) 
"•"'-'-0.01(0.01) 

n n^+o-06(o.32) 

"■"3-0.01(0.01) 
„ ,0+0. 30(^) 
"■^°-0.14(0.14) 



0.8 
0.7 
0.9 
0.7 



Low dust extinction; Ry = 1 



HB06c 

Y08 S(0) 

O081 5(0) 

O08u 5(0) 



Cole et al. 



1 2001 1 with values from HB06 
3.4, zt = 1, 72 = -0.3 

3, Zj, = 1, 72 = -2 

4, zt, = l, 72 = 



17.8, 71 
17.8, 71 
17.8, 71 



1.02 
0.99 



+0.10(0.20) 
-0.09(0.19) 
+0.09(0.18) 
-0.08(0.17) 
, , 0+0.13(0.28) 
■ °-0.22(0.26) 

n Qn+0-08(0-15) 
"•""-0.07(0.15) 



0.8 
0.7 
0.7 
0.7 



"■" -0.00(0.00) 
04-^0-«(l-24) 

"■"^-o.oo(-o.oo) 

n n^+0.06(0.29) 
"■"3-0.01(-0.01) 

+0.21(/) 
"■°'^-0.29(-0.51) 



0.8 

0.7 
0.8 
0.8 



"<!/(<) (X f-i/^ po„er law at f < tc, and ■jfjt) oc a t t > t 
''SFH references: HB06c - 



Yiiksel et al. 



(2008 I ; O081 and O08u - 



Oda et al. 



77i 



I 2008 I 



Hopkins & Beacom (20061; Y08 
Except for HB06c, all other Sl'Hs are parametrized as broken power laws ot the form 5(z) = 5(0) (T 
with 71 at 2 < zi,, and 72 at 2 > 25. 5(0) is in units of 10~^ Mq yr"-' Mpc" ^. 

''The first and second errors (in parentheses) are 68 and 95 per cent confidence regions, respectively, for the slope fH 
of the power-law DTD. 

■^The first and second errors (in parentheses) are 68 and 95 per cent confidence regions, respectively, for tc, the break 
between a t~^/^ and a power law. 

fAs the O08u SFH was found to be compatible with a t power-law DTD at all times, there is no 95 per cent 
upper limit for this measurement. 



mal X? = 1-1 for 35 DOF. However, all of the other SFHs 
result in best-fitting values with minimal Xr > 1-5, ruling 
out this model at the 95 per cent confidence level. The re- 
sulting SN la rate evolution tracks are plotted in the bottom 
panel of Fig. [U] 

We propagate the systematic uncertainty brought about 
by the possibility that the extinction law in the immediate 
vicinity of SNe is different from the Galactic average by fit- 
ting the different DTDs to the SN la measurements derived 
with Ry = 1, as detailed in Section |7.1.2[ The [Pritchet] 
power-law and DOS Gaussian DTDs are 



et al. 



(20081 t 



-1/2 



still excluded, at the 95 per cent confidence level, when using 
the same SFHs as detailed above. The resulting best-fitting 
parameter for the power-law and broken DTDs are pre- 
sented in Table [7] The lower Ry value adds a systematic un- 
certainty of lo'57 to the best-fitting value of /3 for the Y08 
SFH. The overall best-fitting value of /3 for the Y08 SFH 
is thus /? = 1.1 ± 0.1(0.2) (statistical) ± 0.17 (systematic), 
where the statistical errors are the 68 and 95 per cent (in 
parentheses) confidence regions, respectively. The upper 68 
per cent limit on for the broken-power-law DTD rises to 
0.76 Gyr, and the 95 per cent upper limit afforded by the 
Y08, HB06c, and O081 SFHs rises to 1.43 Gyr. 



9 THE TYPE lA SUPERNOVA RATE AT 
REDSHIFT > 2 AND COSMIC IRON 
ACCUMULATION 

Our analysis, above, has provided the most precise deter- 
mination to date of the SN la rate at 1 < z < 2. As seen 
in the bottom panel of Fig. |13[ the best-fitting power-law 



DTD, convolved with each SFH, can also be used to pre- 
dict the SN la rate at z > 2. The shaded regions in the 
figure show the uncertainty regions produced by the statis- 
tical and systematic uncertainties of the DTD slope /3, where 
the statistical uncertainties result from the SN la rate mea- 
surements, and the systematic uncertainties result from the 
uncertainty in the slope of the SFHs at z < zt- 



Following [Blanc fc G rcggio (2008), we can use our re- 
sults to calculate the mean cosmic accumulation of iron. A 
typical SN la produces ~ 0.7 Mq of iron (e.g., Mazzali et al. 
20071. We integrate over the SN la rate evolution derived 



from convolving the power-law DTD described in the previ- 
ous section with the Y08 SFH, multiplied by the above iron 
yield, to derive the amount of iron produced by SNe la. The 
uncertainty in the amount of iron contributed by SNe la is 
calculated by integrating the upper and lower bounds of the 
shaded area in Fig. |13[ multiplied by the above iron yield. 
This takes into account both the spread in the SN la rate 
measurements, and the plausible range of SFH shapes. We 
calculate the amount of iron produced by CC SNe by in- 
tegrating over the Y08 SFH fit. Using the |Salpete7] ( |1955[ ) 
IMF (as assumed by Y08), we calculate either the number 
of stars with masses 8 < M < 50 or the mass in such 
stars. If we assume that 1 per cent of the CC SN progenitor 
mass is converted into iron (as in Maoz et al. 20101, then 



the present-day ratio of iron mass produced by SNe la to 
that produced by CC SNe is 1:4. If, on the other hand, we 
assume that each CC SN produces, on average, 0.066 Mq of 
iron (as in Blanc fc Greggio|2008 based on CC SN samples 
from |Zampieri et al.|2003| and |Hamuy |2003[ ), then the ratio 
increases to 1:1. As the major source of uncertainty in the 
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Figure 15. Cosmic iron density as a function of redshift. In both panels the SN la contribution is denoted by the dot-dashed line, the 
CC SN contribution by the dashed line, and the total amount of iron by the solid line. The dark region around the SN la contribution 
is the systematic and 68 per cent statistical uncertainty introduced by the SFH fits and the SN la rate measurements, respectively. The 
shaded region around the CC SN contribution is the result of the systematic uncertainty in the SFH fits alone. The dark region around 
the total iron density curve is the uncertainty introduced by both SN components. Thin lines delineate the uncertainty regions of each 
component. Left: assuming 1 per cent of the CC SN progenitor mass is converted into iron. Right: assuming each CC SN, on average, 
produces 0.066 Mq of iron. 



amount of iron contributed by CC SNe is the normalization 
of the SFH, we integrate over the O08u and HB06c SFHs 
to derive upper and lower bounds on the uncertainty region. 
Finally, we sum the lower (upper) uncertainty bounds of the 
separate SN la and CC SN contributions to arrive at lower 
(upper) limits on the total cosmic density of iron. 

Both scenarios are presented in Fig. |15[ The mean cos- 
mic iron abundance in solar units, marked on the left ordi- 
nate axis, is calculated assuming i^t — 0.0445 for the baryon 



densi ty in units of the critical closure density ( Komatsu et al. 
201l[ ), and Zpe.p = 1.3 ±0.1 x 1 0"^ for the solar iron abun- 
dance i Grevesse &: Sauval|1998 L We see that the predicted 
present-day mean cosmic iron abundance lies in the range 
(0.09-0.37) ZFe,Q- Between 2 = and 2, for a given choice 
of SFH, the abundance behaves roughly linearly as 



Zfc,L ' 



0.36 - 
! 0.20 



0.10(1-^2), 

-0.06(1 + 2), 



(10) 



for the best-fitting solid black curves in the left (L) and right 
(R) panels of Fig. 15 respectively. The choice of SFH prop- 
agates to a dominant systematic uncertainty in the CC SN 
contribution to the iron abundance. 

From the figures above, we see that the bulk of the pre- 
dicted IGM enrichment occurs at 2 < 2. At these epochs, 
most of the IGM (holding the majority of baryons in the 
Universe) is in the warm-hot intergalactic medium (WHIM) 
phase. The WHIM has yet to be detected clearly in the X-ray 
absorption lines of intermediate elements, let alone of iron, 
which is extremely challenging. However, it is conceivable 
that future X-ray missions, such as the International X-ray 
Observatory ( [Barcons et al.|[2011[ ) or the recently cancelled 
EDGE, having large effective areas and high spectral res- 
olution, could detect FeXVII absorption at A « 17A, and 
eventually lead to an actual low-2 iron abundance measure- 
ment ( Paerels et al.|2008 l. Such a measurement can then be 



compared with these predictions to constrain both the in- 
tegrated iron production of CC SNe and the efficiency with 
which metals produced by SNe are ejected into the IGM. 



10 CONCLUSIONS 

By surveying four deep epochs of the 0.25 deg^ SDF, we 
have assembled a sample of 150 SNe, of which 26 are SNe la 
at 1.0 < 2 < 1.5, and 10 are SNe la at 1.5 < 2 < 2.0. 
This is the largest sample of SNe la at such high redshifts 
to date. The number of SNe la in our 1.0 < 2 < 1.5 bin 
is comparable to that of DOS in the same range, but our 
1.5 < 2 < 2.0 sample is 2.5 times as large. While we may 
have discovered some non-la transients in the redshift range 
1.5 < 2 < 2.0, we have argued that further contamination of 
our high-2 SN la sample is unlikely. Through various tests, 
we have shown that the high-2 SNe in our sample are se- 
curely associated with galaxies at these redshifts, and since 
our survey is mostly insensitive to CC SNe, they must be 
SNe la. The SN la rates derived from our sample are consis- 
tent with those of DOS, but are 2-3 times more precise, with 
uncertainties of 30-50 per cent. Our measurements indicate 
that, following the rise at < 2 < 1, the SN la rate appears 
to level off after 2 ~ 1, but there is no evidence for a decline 
in the SN la rate evolution of the form advocated by DOS. 

Based on these rates and on a growing number of ac- 
curate measurements at 2 < 1, and combined with different 
SFHs, we find that a power-law DTD of the form *(*) = 
*i(t/l Gyr)" fits the data weU, with /3 = -1.1±0.1(0.2) (6S 
and 95 per cent statistical confidence, respectively) ±0.17 
(systematic). This form is consistent with the DTDs found 
by most of the recent SN la surveys, in a variety of envi- 
ronments, at different redshifts, and using different method- 
ologies flbtani et al.|[2008{ |Maoz et al.][20TT| [20To} |Maoz| 



fc Badenesff2010[ ). A t ''^^ power law at all delay times, as 
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proposed by Pritchet et al. (20081, is marginally consistent 



with the data. DTDs consisting of broken power laws are 
also acceptable, as long as tc, the time at which the DTD 
transitions from a t~^^^ power law to a power law, is 
less than ~ 0.8 Gyr (68 per cent confidence). The Gaussian 
DTD proposed by D04, D08, and|Strolger et al.|p004l|2010 1 



is ruled out by all but one of the SFHs tested here. Overall, 
these results are suggestive of the DD progenitor scenario, 
for which a power law with /3 « — 1 is a generic prediction. 
In contrast, DTDs predicted by SD models have a variety 
of forms, but as a rule, they fall off steeply or cut off com- 
pletely beyond delays of a few Gyr (e.g., Meng, Li, & Yang 



2011). The DD channel may not be the only one that pro- 



duces SNe la, but it appears that a large fraction of SNe la 
are formed in this way, or in some other way that mimics 
the DTD predictions of the DD channel. 

Using the best-fitting power-law DTD, we have recon- 
structed how the mean iron abundance of the universe has 
evolved with cosmic time, and predict it is now in the range 
(0.09-0.37) ^Pe,0. This prediction is consistent with those 
of [Fukugita fc Peebles] ^0Q^ and [Blanc fc Greggio] ( [20081 ) , 
but is now based on the most recent and accurate SN la rate 
measurements, the full range of plausible cosmic SFHs, and 
the current DTD estimations. 

The time-integrated number of SNe la per unit mass 
derived from the best-fitting power-law DTD, assuming 
a 'diet-Salpeter' IMF ( [Bell et aL[[2003[ ), is in the range 
Nsn/M^ = (0.5-1.5) X lO""* Mq , though it might eas- 
ily be higher if the normalization of the SFH is found to be 
lower than currently assumed. 

The CC SN rate at (z) = 0.66 is 6.9l^ j x 10"** yr"* 
Mpc~^. This value is consistent with the only other mea- 
surement in this redshift range (D04), and shows that, as 
expected, the CC SN rate tracks the cosmic SFH out to 
2 « 1. 

Our survey in the SDF has reached the point where the 
systematic uncertainties in the SN rates are comparable to 
the statistical uncertainties. The 1.5 deg^ Hyper-Suprime 
Cam (Furusawa et al. 2010), soon to be installed on the 



Subaru Telescope, could allow discovery of larger numbers 
of SNe per epoch and thus a further reduction in the statis- 
tical uncertainties. A new SN survey in a well-studied field, 
such as the SDF or the SXDF, but with cadences designed 
to probe the light curves of the SNe, could permit classifica- 
tion of the SNe at a higher level of accuracy, thus reducing 
the systematic uncertainties as well. This will also apply 
to future massive surveys such as the Large Synoptic Sur- 
vey Telescope ( Stubbs et al.|[2004 ) or the Synoptic All-Sky 
Infrared Survey (Bloom et al. 2009), for which traditional 



spectroscopic foUowup observations will be impossible, but 
to which the approach we have adopted here is optimally 
suited. 

Two //ST Treasury programmes — CLASH (GO12065- 
12069, GO12100-12104) and CANDELS (GO12060-12061) 
— have recently begun deep IR observations utilizing the 
F125W and F160W filters on the Wide Field Camera 3. 
These filters, similar to J and H, will probe the optical part 
of the SN spectrum out to z ~ 1.5, and the near-UV part of 
the spectrum out to 2 ~ 2.7. By observing the optical part of 
the spectrum in the observer-frame IR, one can reduce the 
uncertainties due to high-redshift dust, thus lowering the 
systematic uncertainty of the SN rates in the redshift range 



1.0 < z < 1.5. Ultimately, these programmes will provide 
independent measurements of the SN la rate in the 1.0 < 
z < 2.0 range probed by this work, as well as extend our 
knowledge of the SN la rate evolution out to 2 ~ 2.7. Based 
on the results presented here, as seen in Fig. |13| we predict 
that CLASH (CANDELS) will observe 10-24 (9-19) SNe la 
at 1.0 < 2 < 2.0, and 0-4 (2-7) SNe la at 2.0 < 2 < 2.7. 
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Table 8. 1.5 < z < 2.0 SNe discovered in the SDF. The full table, including the entire sample, is available in the electronic version of the paper. 
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(1) — SN identification. 

(2) -(3) - Right ascensions (J2000; starting at 13*") and declinations (J2000; starting at +27°). 

(4) - SN offset from host galaxy, in arcseconds. Uncertainties appear in parenthesis, and have been multiplied by 100. 

(5) -(7) - SN photometry in the R, i' , and z' bands, in magnitudes. Uncertainties appear in parenthesis, and have been multiplied by 100. 

(8) - Signal-to-noise ratio of the SN, as measured in the z'-band image. 

(9) -(10) - Photometric redshift of SN host galaxy, with reduced , as derived with ZEBRA. 

(11) - Spectroscopic redshift of SN host galaxy, where available. 

(12) -(14) - Probability of a SN being a Type la, or CO SN, as derived with the SNABC, together with its posterior redshift and reduced x^- 
(15)-(16) - Final adopted SN type and redshift. 



Table 9. 1.5 < z < 2.0 SN host galaxies. The full table, including the entire sample, is available in the electronic version of the paper. 
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05(05) 


23.86(06) 


23.94(06) 


24.37(10) 


23.19(18) 


23 


16(14) 


hbUt U7U5 


29 


25 


01 


78 


18 


38 


96 







24 


17(02) 


23 


95(04) 


23 


59(03) 


23 


12(03) 


z2.71(U2j 


22.75(02) 


22.74(02] 




20 


81(03) 


hSDF0806 


31 


24 


19 


54 


29 


59 


51 







26 


17(09) 


25 


85(14) 


25 


42(09) 


24 


69(07) 


24.01(06) 


24.56(10) 


24.11(08) 




20 


80(03) 


hSDF0806 


32 


25 


20 


46 


43 


08 


35 







25 


74(06) 


25 


21(09) 


25 


43(10) 


25 


46(12) 


25.17(15) 


26.62(34) 


> 26.54 








hSDF0806 


38 


23 


33 


35 


14 


20 


69 








24 


93(03) 


24 


29(02) 


24 


23(02) 


23 


93(02) 


23.82(04) 


23.56(03) 


23.65(04) 


22.81(15) 


22 


30(08) 


hSDF0806 


46 


24 


29 


98 


14 


09 


13 







27 


23(17) 


25 


87(14) 


25 


75(12) 


24 


58(07) 


23.83(06) 


24.10(07) 


23.78(06) 




21 


23(04) 


hSDF0806 


50 


23 


46 


02 


38 


59 


59 







25 


36(05) 


24 


56(05) 


24 


28(04) 


23 


95(04) 


23.15(03) 


23.91(06) 


23.23(04) 


21.68(08) 


20 


61(03) 


hSDFOSOe 


57 


25 


33 


67 


28 


03 


27 







25 


66(06) 


25 


65(12) 


25 


16(08) 


24 


88(08) 


24.47(09) 


24.78(12) 


24.77(13) 




22 


51(08) 



Note - magnitude limits are 3<T. 

(1) - SN identification. 

(2) -(3) - Right ascensions (J2000; starting at IS*") and declinations (J2000; starting at +27°). 

(4)-(5) - GALEX FUV axid Afi/y photometry. —1 means no UV signal observed in this bajid; 1 means a clear UV signal associated with the target galaxy; 

and means the UV signal could not be unequivocally matched to the target galaxy. 
(6)-(12) - Subaru optical photometry, in magnitudes. Uncertainties appear in parenthesis, and have been multiplied by 100. 
(13)-(14) - UKIRT J and K photometry, in magnitudes. Uncertainties appear in parenthesis, and have been multiplied by 100. 



